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ABSTRACT 


A  technique  is  presented  for  predicting  the  vibration  responses  to 
reentry  aero-acoustic  loading  of  high  beta  ballistic  vehicles  (those 
having  a  weight-to-drag  ratio  greater  than  1500).  Correlation  of 
predicted  responses  due  to  ground  acoustic  loading  and  measured  data 
is  given  and,  based  on  these  results,  guidelines  on  the  use  of  the 
prediction  technique  for  acoustic  and  aero-acoustic  loading  are 
outlined.  Computer  programs  for  calculating  vibration  responses  due  to 
reverberant  and  progressive  wave  acoustic  environments  and  to  aero- 
acoustic  excitations  are  documented  and  arc  available  through  the 
\i r  Force  Flight  Dynamics  Laboratory  (FYj.  Techniques  for  ground 
simulation  of  re-entry  vibration  response  are  discussed,  and  typical 
component  design  requirements  and  test  specifications  are  delineated. 
The  new  feature  of  the  technique  Is  the  use  of  finite  element  methods 
in  deriving  component  test  specifications.  These  more  accurate  test 
specifications  will  he  essential  in  reducing  the  number  of  overdesigned 
components,  thereby  decreasing  the  cost  of  mission  completion. 


TABLE  OF  CONTENTS 


Section  Title  Page 


I  INTRODUCTION  . .  1 

II  SUMMARY .  2 

1.  Typical  Re-Entry  Environment .  2 

2.  Response  Prediction  Technique .  2 

3.  Prototype  Vehicle  Dynamic  Characteristics .  2 

4.  Re-Entry  Vibration  Response  Predictions .  4 

5.  Response  Verification  Test .  4 

6.  Prediction  Technique  Verification .  4 

7.  Simulation  Techniques,  Design  and  Test  Requirements .  5 

HI  TYPICAL  RE-ENTRY  ENVIRONMENT .  6 

1.  Aeroacoustic  Environments  .  6 

2.  Forebody  Pressure  Cross  Power  Spectral  Density .  6 

3.  Wake  Pressure  Cross  Power  Spectral  Density .  9 

4.  Octave  Band  Sound  Pressure  Levels .  13 

IV  RESPONSE  PREDICTION  TECHNIQUE .  14 

1.  Literature  Survey .  14 

2.  Optimized  Technique .  16 

3.  Re-Entry  Vibration  Response .  22 

4.  Reverberant  Chamber  Response .  27 

5.  Progressive  Wave  Chamber  Response .  28 

6.  Computer  Programs  .  29 

V  PROTOTYPE  VEHICLE  DYNAMIC  CHARACTERISTICS . 31 

1.  Definition  of  Prototype  Vehicle  .  31 

2.  Dynamic  Model .  31 

3.  Modal  Survey  Summary  .  34 

4.  Frequencies,  Mode  Shapes  and  Damping .  40 

VI  RE-ENTRY  VIBRATION  RESPONSE  PREDICTIONS . 47 

1.  Combined  Forebody  and  Wake  Environments  .  47 

2.  Separate  Forebody,  Wake  Environments .  55 


preceding 


page  blank-not  fumed 


V 


TABLE  OF  CONTENTS  (Continued) 


Section  Title  Page 


VII  RESPONSE  VERIFICATION  TEST .  56 

1 .  Objcelh  cs  .  56 

2.  Test  Requirements .  5(5 

3.  rest  Results . 61 

VI fl  PREDICTION  TECHNIQUE  VERIFICATION .  6(5 

1 .  Test  Acoustic  Environment  Description .  66 

2.  Test  Vibration  Response  Prediction .  72 

6.  Test- Analysis  Correlation .  72 

I.  Discussion .  79 

IX  VEHICLE  RE-ENTRY  SIMULATION  TECHNIQUES .  “0 

1.  Current  Aeroacoustic  Environment  Description .  >0 

2.  Current  Re-Entry  Vibration  Environment  Description .  32 

2 .  Simulation  Techniques .  72 

4.  AFFDL  Sonic  Fatique  Facility  Capabilities .  «5 

X  DESIGN  REQUIREMENTS,  TEST  SPECIFICATIONS  .  3(5 

1 .  Specification  Derivation  Guidelines .  SO 

2.  Typical  Component  Test  Specifications .  37 

3.  Internal  Acoustic  Envii’onments . . .  37 

XI  CONCLUSIONS  AND  RECOMMENDATIONS  .  79 

1 .  Vibration  Response  Predictions .  39 

2.  Acoustic  Testing .  69 

3.  Simulation  Techniques  .  90 

REFERENCES  .  91 


Figure 


&> 

fcj 

ft  * 


b 


1 

2 

3 

4 

5 


i 


i 


6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 


LIST  OF  ILLUSTRATIONS 


Title 


Summary  of  Vibration  Response  Prediction . . . 

Flow  Parameters  Along  the  Prototype  Vehicle  . 

Spatial  Correlation  Functions . 

Fluctuating  Aero-Dynamic  Properties  on  Base  of  Re-entry  Vehicle  .  .  .  . 

Delta  Turning  Angle  Vs,  Mach  Number . 

Octave  Band  Soiuid  Pressure  Levels  on  Forebody  and  Base . 

SABOR  Finite  Element  Coordinate  Conventions . 

Dynamic  Model . 

Sabor  Modeling  of  Shell  Joints . 

Sabor  Modeling  of  Bulkhead  Attachments . 

Modal  Survey  Test  Setup . 

Schematic  of  Sensor  Identification  and  Location . 

Axial  Mode  Shape . 

Axial  Mode  Shape . 

Lateral  Me-***  Shape . 

Lateral  Moo-,1  Shape  . 

Lateral  Mode  Shape . 

Accelerometer  1AI,  Response  . 

Accelerometer  1AT  Response  . 

Accelerometer  3AL  Response  . 

Accelerometer  3 AT  Response  . 

Accelerometer  4AL  Response  . 

Accelerometer  4 AT  Response  . 

Accelerometer  5AL  Response  . 

Accelerometer  5AT  Response  . 

Accelerometer  6AL  Response  . 

Accelerometer  6AT  Response  . 

Accelerometer  8AN  Response  . 

Accelerometer  9AN  Response  .  . . 

Accelerometer  10 AL  Response . 

Vehicle  Support  System  . 

Schematic  of  External  Acoustic  Sensors . 

Acoustic  Spectrum  —  Minimum  Level . 

Quick- Look  Acoustic  Data  Processing  . 

Detailed  Acoustic  Data  Pi’ocessing  . 

Vibration  Data  XTocessing . 

Forebody  Sound  Pressure  Level  in  Progressive  Test . 

Forebody  Sound  Pressure  Level  in  Reverberant  Test . 

Base  Sound  Pressure  Level  in  Progressive  Test  . 

Base  Sound  Pressure  Level  in  Reverberant  Test  . 

Internal  Sound  Pressure  Level  (2PL)  in  Progressive  Test  . 

Internal  Sound  Pressure  Level  (7PL)  in  Progressive  Test  . 


vii 


LIST  OF  ILLUSTRATIONS  (Continued) 


Figure  Title  Page 

43  Narrow  Band  PSD'S  of  Microphone  LI .  67 

44  Meridional  Space-Time  Correlation  Data .  67 

45  Meridional  Spatial  Correlation  Function .  68 

46  Circumferential  Space-Time  Correlation  Data .  70 

47  Circumferential  Spatial  Correlation  Function  .  70 

48  Meridional  Space-Time  Correlation  Comparison .  71 

49  Circumferential  Space-Time  Correlation  Comparison .  71 

50  1 AL  Vibration  Response,  Test  Vs.  Analysis  .  73 

51  lAT  Vibration  Response,  Test  Vs.  Analysis  .  73 

52  3 AL  Vibration  Response,  Test  Vs.  Analysis  .  73 

53  3 AT  Vibration  Response,  Test  Vs.  Analysis  .  73 

54  4 AL  Vibration  Response,  Test  Vs.  Analysis  .  74 

55  4 AT  Vibration  Response,  Test  Vs.  Analysis  .  74 

56  5 AL  Vibration  Response,  Test  Vs.  Analysis  .  74 

57  SAT  Vibration  Response,  Test  Vs.  Analysis  .  74 

58  6 AL  Vibration  Response,  Test  Vs.  Analysis  .  75 

59  6 AT  Vibration  Response,  TescVs.  Analysis  .  75 

60  8 AN  Vibration  Response,  Test  Vs.  Analysis  .  75 

61  9 AN  Vibration  Response,  Test  Vs.  Analysis  .  75 

62  10AL  Vibration  Response,  Test  Vs.  Analysis  .  76 

63  .Axial  Diode  Shape,  Test  Vs.  Analysis  .  78 

64  Current  Aeroacoustic  Environment  Description .  81 

65  PSD  Envelope  for  Internal  Components-Transition .  82 

66  PSD  Envelope  for  Internal  Components  -  TBL  .  83 

67  Re-entry  Response  Simulation  Technique  .  84 

68  Test  Specification  Notching .  86 

69  Component  Re-entry  Vibration  Test  specification  (Internal  Self 

Mounted)  .  87 

70  Component  Re-entry  Vibration  Test  Specification  (Forebodv  Shell 

and  Base)  .  88 

71  Component  Re-entry  Shock  Test  Specification  .  88 


LIST  OF  TABLES 


Table  Title 

I  Coordinates  at  Nodes . 

II  Weight  Breakdown . 

III  Sensor  Description  . 

IV  Comparison  of  Axial  Experimental  and  Analytical  Frequencies  . 

V  Comparison  of  Lateral  Experimental  and  Analytical  Frequencies 

VI  Transformations  for  Accelerometer  Response  Predictions  .  .  .  . 


Page 

33 

34 
39 

41 

42 
54 


LIST  OF  SYMBOLS 


A„  —  Area  of  n^1  element,  ft^ 

a*  —  Speed  of  sound  where  Much  number  equals  1,  ft/sec 

C  —  Co  narrowband  space  correlation  coefficient 

c  —  Speed  of  sound  in  sonic  chamber,  ft/sec 

Ci,  c2  —  Decay  constants  in  co-space  correlation  coefficient 

f  —  Frequency,  IIz 

fc  —  Center  frequency,  Hz 

fj-  —  Frequency  bandwidth,  Hz 

g;  —  Structural  damping  coefficient 

H|.  —  Modal  transfer  function,  with  +  denotes  complex  conjugate 

i  —  Imaginary  designation 

J  ,  q  —  Joint  acceptance  in  circumferential  direction 

K  —  Stiffness  in  typical  equation  of  motion 

i  L  1  ~  Displacement  transformation  matrix 

?  —  Meridional  length  of  finite  element,  ft 

M  —  Mass  in  typical  equation  of  motion,  Ib-sec^/ft 

M^  —  Mach  number  associated  with  base  flow 

Mc  —  Local  or  edge  Mach  number 

M|c  —  Modal  mass  of  mode,  lb-sec“/ft 

n  —  Circumferential  wave  number 

9 

p  (s,  t)  —  Distributed  pressure  in  typical  equaHon  of  motion,  Ib/ft“ 
P^  —  Base  static  pressure,  lb/ft2 


P^SWSSsi«(3S--s9(t 


PT 

{SI 
If  I 

Q 


c4 


R 


sT 
!  ) 


V, 


w 


|i 

iP 


'  f*#*6"  ^v^T?V 


Total  pressure,  lb/l't- 

Generalized  forces  on  finite  element  in  element  coordinate  system 
Generalized  forces  on  finite  element  in  system  coordinates 
Quadrature  narrowbard  space  correlation  coefficient 
Dynamic  pressui-e,  lbs/ft“ 

Finite  element  coordinates 
Radius  of  base,  ft 

Space -time  cross  correlation  function 
Radius,  ft 

Radius  on  number  1  finite  element,  ft 
Meridional  location,  ft 

Meridional  distance  on  number  1  finite  element,  ft 

Meridional  separation  lx?  tween  left  station  of  elements  1  and  2,  ft 

Transpose  of  matrix  in  brackets 

SABOR  axial  displacement 

Convected  velocity  of  turbulent  flow,  ft  ./sec 

Local  or  edge  flow  velocity,  ft/sec 

SABOR  tangential  displacement 

SABOR  radial  displacement 

Coordinate  in  typical  equation  of  motion 

Vehicle  cone  half  angle,  degrees 

SABOR  edge  rotation 

Transformation  matrix  to  convex!  finite  elemental  coordinates  to  system 
coordinates 


1 


i 


1  \» 
J  3 


r! 


* 


> 


'  i 


Modal  damping 


Damping  in  typical  equation  of  motion 
Boundary  layer  displacement  thickness,  ft 
Separation  distance  across  flow  dix-ection,  ft 
Circumferential  direction  location,  longitude  angle,  rad 
Half  angle  of  wake  cone,  degrees 
Circumferential  direction  location  on  element  1,  rad 
Prandtl-.Mayer  expansion  angle  on  base,  degrees 
Local  Prandtl-Mayer  expansion  angle,  degrees 
Separation  distance  along  flow  direction,  ft 
Modal  coordinate  in  typical  equation  of  motion 
Generalized  force  integral  expressions 

3 

Local  or  edge  density  lb/ft 
Base  rms  pressure,  lb/ft2 
Local  rms  pressure,  lb/ft2 
Delay  time,  sec 

Base  pressure  cross  power  spectral  density,  psfVlIz 
Pressure  cross  power  spectral  density  for  turbulent  flow,  psfvHz 

ft 

Pressure  cross  power  spectral  density,  psf  /Hz 

Force  cross  power  spectral  density  lb2/Hz 

Co-force  cross  power  spectral  density,  lbVflz 

Quad  -  force  cross  power  spectral  density’,  lb2/ilz 

Term  in  force  cross  power  spectral  density  matrix,  lbVliz 


Response  cross  power  spectral  density,  response  unitsVliz 

Displacement  response  power  spectral  density,  j  2/IIz 

Homogeneous  pressure  power  spectral  density.  Subscripts  denote 
location,  psfVlIz 

Mode  shape 

Incidence  angle  for  progressive  waves,  degrees 
Radian  frequency 

Radian  frequency  of  natural  modes  of  vibration 


SECTION  I 

INTRODUCTION 

The  prediction  of  vibration  responses  induced  by  random  pressure  fluctuations  during  re¬ 
entry  of  a  high  beta  vehicle  is  a  difficult  task.  (Beta  is  the  ratio  of  vehicle  weight  to  the 
product  of  its  drag  coefficient  times  its  maximum  cross-sectional  area. )  In  general,  the 
work  of  industry  has  been  predominantly  theoretical  with  very  limited  experimental  data 
available  to  adequately  determine  the  reliability  of  these  theoretical  methods.  As  a  result  of 
low  confidence  in  predictions,  it  has  been  necessary  to  impose  very  conservati  e  vibration 
design  and  test  criteria  in  the  development  of  components  for  high  beta  re-entry  vehicles. 
These  test  specifications  are  usually  obtained  by  multiplying  scant  measured  flight  data  by 
the  ratio  of  the  dynamic  pressures  between  the  flight  vehicle  and  the  design  vehicle  for  which 
the  specifications  are  intended.  This  technique  omits  many  factors  which  influence  the 
vibi’ation  response  and  thus  results  in  excessively  high  specifications.  Thus  much  contract 
time  and  funding  is  expended  in  overdesigning  components  to  pass  laboratory  tests  that  bear 
little  resemblance  to  actual  re-entry  environments. 

The  objective  of  this  exploratory  and  development  program,  then,  is  to  assess  the  accuracy 
of  analytical  methods  used  throughout  industry,  and  to  refine  and  verify  a  selected  method  to 
enable  the  establishment  of  optimum  vibration  design  and  test  critex-ia  for  both  structure  and 
equipment.  These  optimum  design  and  test  criteria  will  be  much  more  representative  of 
actual  flight  environments  thereby  reducing  component  specification  oi*r  conservatism. 
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SECTION  II 


SUMMARY 


The  prediction  of  re-entry  vibration  response  involves  the  use  of  aeroacoustic  forcing  func¬ 
tions  in  conjunction  with  vehicle  dynamic  characteristics.  The  description  and  interaction 
of  these  input  parameters  with  the  prediction  technique  is  illustrated  by  the  flow  diagram  in 
Figure  1.  This  flow  diagram  forms  the  basis  for  the  organization  and  logic  sequence  of 
tasks  performed  in  this  study  program;  in  particular,  it  emphasizes  the  need  for  well- 
defined  acoustic  forcing  functions  and  vehicle  dynamic  characteristics  in  order  to  achieve 
verification  of  the  response  prediction  technique. 

1.  TYPICAL  RE-ENTRY  ENVIRONMENT 

Utilizing  state-of-art  techniques  in  defining  fluctuating  pressure  environments  on  re-entry 
vehicles,  a  baseline  re-entry  environment  is  defined  for  the  prototype  vehicle  to  be  used  in 
later  tasks.  There  are  two  ae rodynamic  induced  environments  considered.  These  are 
pressure  fluctuations  in  the  turbulent  boundary  layer  which  act  on  the  vehicle  forebody,  and 
pressure  fluctuations  in  the  wake  which  act  on  the  vehicle  base.  The  homogeneous  pressure 
spectrum  is  defined  at  each  vehicle  station  as  well  as  cross  correlation  coefficients  so  that 
the  cross  pressure  spectrum  can  be  represented.  The  aerodynamic  parameters  utilized  in 
these  expressions  included  local  Mach  numbers,  local  convected  velocities,  local  densities, 
and  boundary  layer  displacement  thickness.  These  aerodynamic  parameters  are  not 
typically  generated  by  a  structural  dynamicist.  Reference  (23)  shows  how  these  parameters 
are  obtained  from  the  information  which  is  usually  available  such  as  trajectory  properties. 

2.  RESPONSE  PREDICTION  TECHNIQUE 

An  all-inclusive  formulation  of  vehicle  response  due  to-distributed  random  pressure  environ¬ 
ments  is  first,  presented.  In  this  formulation  no  simplifying  assumptions  have  been  made. 
Then,  utilizing  the  results  of  a  detailed  literature  survey,  the  complex  formulation  is 
simplified  to  yield  a  less  complex  response  prediction  technique,  with  no  appreciable  loss  in 
accuracy  but  significant  gains  in  economy  associated  with  the  computer  implementation.  A 
detailed  examination  of  the  force  cross  Power  Spectral  Density  matrix  is  made  to  show  how 
the  analytical  techniques  can  be  applied  to  acoustic  environments  encountered  in  ground 
testing  as  well  as  re-entry  aeroacoustic  environments.  It  is  shown  how  further  simplifica¬ 
tions  can  be  made  to  the  formulation  of  the  force  cross  PSD  matrix  in  cases  where  the 
pressure  cross  PSD  can  be  considered  constant  over  the  size  of  the  finite  elements  in  the 
dynamic  model.  In  the  more  general  case  where  this  assumption  is  not  valid,  individual 
terms  in  the  force  cross  PSD  matrix  are  shown  as  they  would  be  numerically  integrated. 

3.  PROTOTYPE  VEHICLE  DYNAMIC  CHARACTERISTICS 

A  mathematical  representation  of  the  dynamic  characteristics  of  the  prototype  vehicle  is 
formulated.  From  available  drawings  and  measurements  of  component  weights  and 
locations,  a  finite  element  idealization  of  the  vehicle  is  obtained.  The  shell  modeling  is 
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ACOUSTIC  INPUT  MODEL 


Figure  1.  Summary  of  Vibration  Response  Prediction 


performed  utilizing  the  MIT  -  developed  SABOR  program  and  the  vehicle  nose  and  component 
packages  are  modeled  using  beam  finite  elements.  Two  separate  dynamic  models  are  con¬ 
structed,  one  describing  the  vehicle  axial  dynamic  characteristics  and  the  other  describing 
the  lateral  dynamic  characteristics . 

Io  verify  and  supplement  the  dynamic  models,  a  modal  survey  is  conducted  on  the  prototype 
vehicle.  The  vehicle  natural  frequencies  and  mode  shapes  are  determined  in  the  20  to 
2000  Hz  frequency  range.  This  modal  data  is  then  used  to  l'evise  the  analytical  model  in 
areas  of  obvious  deficiency.  The  modal  s.urvey  data  is  also  utilized  to  obtain  the  damping 
characteristics  in  each  mode. 

In  addition  to  the  external  probe  accelerometer  used  to  define  the  shell  motion,  thirteen 
internal  accelerometers  are  also  utilized.  These  accelerometers  are  located  in  both  compo¬ 
nent  packages,  the  internal  shell  structure  and  the  aft  cover.  It  is  this  internal  instrumenta¬ 
tion  that  is  used  in  the  acoustic  test  to  obtain  the  component  responses  and  verify  the 
response  prediction  techniques . 

4.  RE-ENTRY  VIBRATION  RESPONSE  PREDICTIONS 

The  response  prediction  technique,  in  conjunction  with  the  verified  dynamic  model,  is  used 
to  predict  the  prototype  vehicle  vibration  response  to  the  baseline  re-entry  environment. 
Acceleratiompower  spectral  densities  are  predicted  at  each  of  thednternal  accelerometer 
locations.  The  vibration  response  predictions  utilized  the  simultaneous  application  of  the 
turbulent  boundary  layer  environment  and  the  wake  environment.  A  discussion  is  also 
included  on  the  effects  of  each  environment  applied  separately.  These  vibration  responses 
form  the  basis  for  acoustic  test  requirements  in  verifying  the  prediction  technique. 

5.  RESPONSE  VERIFICATION  TEST 

Acoustic  tests  are  conducted  at  the  AFFDL  Sonic  Fatigue  Facility  on  the  prototype  vehicle. 
The  objective  of  these  tests  is:  (1)  to  verify  the  response  prediction  technique  by  obtaining 
vehicle  response  data  in  a  form  suitable  for  correlation  with  analysis  at  levels  similar  to 
re-entry,  and  (2)  to  define  vehicle  simulation  techniques  by  making  accurate  measurements 
of  external  and  internal  acoustic  levels  to  define  and  express  experimental  forcing  functions. 

Response  predictions  were  made  on  the  prototype  vehicle  under  ideal  reverberant  and 
progressive  wave  conditions,  which,  when-correlated  with  the  re-enti'y  vibration  predictions, 
defined  a  minimum  sound  pressure  level  requirement  for  the  acoustic  tests  to  simulate  re¬ 
entry  as  best  as  possible. 

6.  PREDICTION  TECHNIQUE  VERIFICATION 

The  measured  acoustic  environments  achieved  in  the  AFFDL  Sonic  Fatigue  Facility  are  used 
in  the  prediction  technique  to  compute  vibration  response  at  the  instrumentation  locations. 
Correlation  of  predicted  responses  with  test  data  shows  good  comparison.  Guidelines  are 
given  in  utilization  of  the  prediction  technique;  in  particular,  attention  is  given  to  dynamic 
modeling  procedures  required  as  a  function  of  the  type  and  location  of  vibration  response 
desired. 
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SECTION  III 


TYPICAL  RE-ENTRY  ENVIRONMENT 


1.  AERO  ACOUSTIC  ENVIRONMENTS 

During  re-entry  a  high  Lota  vehicle  is  subjected  to  pressure  fluctuations  in  the  attached  tur¬ 
bulent  boundary  layer,  pressure  fluctuations  in  the  wake  flow,  and  pressure  fluctuations  due 
to  the  transition  from  laminar  to  turbulent  flow.  If  the  vehicle  is  of  the  class  containing  pro¬ 
tuberances  or  flare  sections,  local  oscillating  shocks  will  also  occur.  The  prototype  vehicle 
used  in  this  study  did  not  contain  protuberances  and  is  therefore  not  expected  to  experience 
shock  loadings  of  any  significance  during  re-entry.  At  this  time,  little  is  known  about  the 
nature  of  the  environment  during  transition  and  there  are  no  analytical  expressions  for 
pressure  cross  PSD  during  transition.  The  environments  presented  for  this  study  are  there¬ 
fore  for  fully  -  turbulent  and  waive  flow  conditions. 

2.  FORE  BODY  PRESSURE  CROSS  POWER  SPECTRAL  DENSITY 

Cl  2) 

Koubolt's  relations  ’  arc  used  to  express  the  pressure  cross  PSD.  There  are  four  flow 
parameters  required  to  formulate  the  pressure  spectrum  and  spatial  correlation  functions. 
These  parameters  are  buundaty  layer  displacement  thickness,  6  *,  local  or  edge  Mach 
number,  Me,  local  or  edge  velocity,  Ve,  and  local  or  edge  density,  pc  .  Figure  2  shows 
how  these  four  parameters  vary  with  prototype  vehicle  station. 

The  pressure  cross  PSD  can  be  written  as 

<f> TBL  ( S ,  if)  =  [ic  ( ij ,  V ,  f)  +  i  Q  M  ,  * .  f)  ]:=*  (f)  (1) 

where 

C(f ,  t?,.f)  -  In-phase  or  co -correlation  function. 

Q(?,tj,I)  -  Out  of  phase  or  quad  correlation  function. 

§  (f)  =  Pressure  PSD  at  a  station. 

X  ~  Separation  distance  in  the  meridional  direction  (along  the 
flow) . 

-  Separation  distance  in  the  circumferential  direction  (across 
the  flow) 
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BOUNDARY  LAYER  DISPLACEMENT  THICKNESS  LOCAL  DENSITY 


Figure  2.  Flow  Parameters  Along  the  Protore  Vehicle 


The  spatial  correlation  functions  can  be  written  as 


where 


“Cll*l  “C2  \  v 

8  *  8  * 
C  ( £ ,  t]  >  f )  =  e  e 


cos 


2tt£{ 

V 

c 


Q  ( $  *  > 


-c 


III 


f)  =e 


sin 


2  ?rf  { 

V 

c 


V  =  Convected  velocity,  equal  to  0.7  Ve. 
c 

c  =  Decay  constant,  assumed  to  be  0 . 05  for  this  study  (based  on  data 
in  .Reference  3) . 

c  =  Another  decay  constant  assumed  equal  to  1^. 

A 

Figure  3  shows  a  graphical  description  of  the  above  correlation  functions . 


(2) 

(3) 


Figure  3.  Spatial  Correlation  Functions 


Many  investigators  also  include  additional  decay  terms  in  the  correlation  functions,  which 
have  the  following  form 


e 


V 

c 


o 


V 

c 
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These  investigators  are  usually  concerned  with  flow  conditions  at  subsonic  or  slightly  super¬ 
sonic  velocities.  For  the  high  beta  re-entry  vehicle  where  convected  velocities  are  10  to  14 
left. /sec  and  the  boundary  layer  displacement  thickness  arc  0.005  to  0.035  ft.,  the  above 


terms  become  insignificant  when  compared  to  '""‘g't"  ~~TT~‘  ^or  examPle>  "the 

total  expression  with  the  added  decay  terms  included  would  be  as  follows:  C(£,  v  ,i)  = 

LI  f  A  J.  -n  I  «  I  J-  "|  2_wf^  TTr.  I  v.  1  ..ntnnr.  £  * 


1 1  r  i  2n-f  i  i  i  r i  2?i-f  I 
v^J  +^JCi 


cos  — —  —  .  Using  stated  values  of  5  * 
'c 


and  V  ,  the  terms  in  the  brackets  would  be 
c 

+4.5  x  10"4fl  . 


1  2?f 

~035  +  14000 


?f 

>00  J.’ 


or  approximately 


It  is  seen  that  for  frequencies  less  than  5  kHz,  the  term  is  negligible.  It  is  for  this  reason 
that  decay  terms  which  are  frequency-dependent  need  not  be  considered  for  structural  re¬ 
sponses  in  the  0  to  5  kHz  range. 

The  pressure  PSD,  ft  (f)  has  been  formulated  by  Houbolt  as: 


*<f)  = 


I  •(¥)■] 


where  a ^  is  the  rms  pressure  equal  to 

a  =  0.00011  P  V  V (1  + 0.012  M  2)  ( 

e  e  e  /  e  ' 

3.  WAKE  PRESSURE  CROSS  POWER  SPECTRAL  DENSITY 

The  analytical  expressions  proposed  for  the  pressure  cross  PSD  over  the  aft  end  of  the 
vehicle  are  essentially  those  for  fully  turbulent  flow  but  with  the  following  modifications: 

(1)  the  boundary  layer  displacement  thickness  is  assumed  equal  to  the  base  radius, 

(2)  the  convected  velocity  is  the  transverse  component  of  the  flow  velocity  along  the 
wake  boundary,  and 

(3)  the  rms  pressure  is  expressed  with  different  coefficients  to  account  for  strong 
viscous  shear  and  lack  of  cooling. 


With  these  modifications  the  following  relations  result: 

‘l>b(i>4,  if)  =  jc($,rj,f)  +  i  Q ( S,  n,  f)]Vf) 
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(?) 

(8) 

(9) 

(10) 


M,  =  Base  Mach  number 
b 

P^  =  Static  base  pressure 

The  effect  of  replacing  the  boundary  layer  displacement  thickness  with  the  base  radius  yields 
a  slowly  decaying  spatial  correlation  function,  indicative  of  highly  correlated  pressures . 
Also  the  R/Vb  sin0*  term  tends  to  cause  the  base  frequency  spectrum  to  decay  much  faster 
than  the  forebody  spectrum. 

The  methods  used  to  determine  the  base  flow  parameters- will  now  be  briefly  described. 

Data  required  to  determine  the  flow  parameters  on  the  base  are  aerodynamic  conditions  on 
the  aft  end  of  the  forebody.  For  this  study  these  properties  are: 

V  =  20325  ft/sec. 

e 

P  =  1759.7  lbs/ft.2 

e 

M  =  6.723 

e 


10 


Tables  giving  properties  of  compressible  flow  including  Prandtl-Meyei  angles,  v  ,  as  a  func¬ 
tion  of  Mach  number  are  required.  An  excellent  source  of  this  information  is  found  m  Ref¬ 
erence  4.  Figure  4  shows  the  flow  parameters  ou  the  re-entry  vehicle  base. 


.•> 


Figure  4.  Fluctuating  Aero-Dynamic  Properties  on  Base 
of  Re-entry  Vehicle 

From  the  tables  of  Reference  4  at  M  =  6.723  we  obtain 

e 

V  =  89.448 
e 

P  /P  =  0.0003127 
e  T 

V  /a*  =  2.3242 
e 

Now  using  Figure  5,  which  is  based  upon  flight  base  pressure  data  ’  and  the  difference  be¬ 
tween  the  Prandtl-Meyer  angles  for  isentropic  expansion  from  cone  pressure  and  base  pre- 
sures^,  we  obtain  A  v  as  16c  at  Mg  =  6.723.  The  Prandtl-Meyer  angle,  associated  with  the 
base  flow  is: 

D  =  V  +  A  v 

b  e 

»  .  =  89.448  4  1(5 
b 

v.  =  105.448 
b 

and  from  Figure  4: 

Q*  -  A  v  -a  =  9.7® 
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TURNING  ANGLE  IN  DEGREES 


LOCAL  MACH  NUMBER 


Figure  5.  Delta  Turning  Angle  Vj.  Mach  Number 

Using  the  tables  of  Reference  4  the  Mach  number,  Mb,  which  eorrespondsXto  a  of 

105.448  is  obtained.  At  this  Mb  the  required  ratios  Pb/PT  vb^a*  are  also  Stained 

M  =11.3 
b 

P/P^  =  0.00001039 
b  i 

V. /a*  =2.403 

b 

We  now  have  all  the  data  required  to  obtain  the  base  flow  parameters. 

p  =  58.54b/ft2 

b  PVPT  e 

V./a* 

Vu  =  ~TT  V  =21014  ft/sec 
b  V  /a*  e 

e 


0.01M  p 

.  - -  =  3.1  lb/ft  rms 

1  +  0.18  M. 


4.  OCTAVE  BAND  SOUND  PRESSURE  LEVELS 


Sound  pressure  levels  for  both  the  base  pressure  spectrum  and  the  maximum  iorebody  spec¬ 
trum  have  been  evaluated  Jn  octave  bands.  Figure  6  shows  the  various  dB  levels  associated 
with  each  octave  band  which  are  used  to  estimate  the  levels  required  in  the  acoustic  test. 
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Figure  6.  Octave  Band  Sound  Pressure  Levels  on  Forebody  and  Base 


SECTION  IV 


RESPONSE  PREDICTION  TECHNIQUE 

I  .  LITERATURE  SURVEY 

The  goal  of  this  task  was  to  obtain,  from  previously  published  work,  simplifying  assumptions 
which  can  be  made  in  the  vibration  response  prediction  technique  with  little  loss  of  accuracy 
but  significant  economic  gains.  Many  reports  and  papers  were  therefore  reviewed  to  satisfy 
this  goal.  Of  particular  interest  in  these  reports  were  the  generalized  force  formulations 
and  the  retention  or  elimination  of  cross  modaL  terms  in  the  response  formulation.  Because 
this  study  made  use  of  a  dynamic  model  constructed  via  the  finite  element  technique,  only 
those  reports  which  also  employed  finite  element  dynamic  models  were  of  interest  with  re¬ 
spect  to  the  gener;  lived  force.  Those  reports  utilizing  closed-form  solution  dynamic  models 
were  of  interest  with  respect  to  the  response  formulation. 

References  6  through  15  are  concerned  with  the  prediction  of  vibration  response  of  simply-^ 
supported  panels  ard  cylinders  to  turbulent  boundary  layer  pressure  fluctuations,  jet  noise, 
and  acoustic  fields.  All  of  the  reports  utilize  a  closed-form  solution  for  the  normal  modes 
of  vibration.  The  response  formulations  nearly  all  follow  the  same  format.  An  all- 
comprehensive  expression  is  initially  developed  which  includes  all  the  cross  modal  contri¬ 
butions  and  the  imaginary’  spectra.  The  formulation  is  then  simplified  to  eliminate  cross 
modal  terms.  There  are  varying  arguments  given  for  this  simplification.  In  the  work  of 
Bozich  and  White  (7'  and  Gcckbcrn  and  Jolly  (9),  space  averaging  the  response  exactly  elimi¬ 
nates  the  cross  modal  response  because  of  orthogonality  of  the  mode  shapes,  Crocker(lO) 
and  Bozich d')  merely  mention  they  are  neglecting  cross  modal  response  terms  and  formulate 
the  response  accordingly,  Maeslrello(H)  (12)  neglects  cross  modal  terms  based  on  Hie 
assumption  that  the  natural  frequencies  are  well  spaced  and -there  is  little  damping  in  the  panel, 

Clarkson(k)  cites  an  investigation  by  another  investigator  as  the  basis  for  elimination  of 
cross  modal  terms  in  his  work.  Tins  other  investigator  showed  only  a  5  to  10  pex*cent  con¬ 
tribution  to  the  rms  response  as  a  result  of  these  terms  in  typical  cases.  SzechenyiD-4) 
gives  as  the  basis  for  elimination  of  cross  modal  terms  the  work  of  Wilby  (15).  Wilby  has 
shown  that  the  cross  modal  terms  do  not  change  the  peak  power  spectral  density  response 
points  at  all.  Wilby  also  shows  good  correlation  with  test  results  in  his  work,  Strawderman 
and  Brand (13)  are  the  exceptions;  they  formulated  the  response  retaining  all  terms. 


It  should  be  pointed  out  that  the  response  formulation  (elimination  of  cross  modal  terms)  was 
the  only  item  applicable  to  this  study  in  References  6  through  15.  However  there  is  much 
more  to  be  gained  from  the  above  reports  if  such  information  as  sound  radiation  from  panels 
and  design  charts  for  fatigue  is  desired. 

References  16  through  19  represent  a  series  of  reports  in  which  the  panel  response  to  turbu¬ 
lent  boundary  layer  pressure  fluctuations  is  obtained  utilizing  a  finite  element  description  of 
the  panel  to  obtain  frequencies  and  modes.  References  16,  17,  and  18  use  the  same  finite 
element  program  and  response  program  for  vibration  solutions.  The  response  formulation 
is  developed  such  that  all  terms  are  utilized  and  both  real  and  imaginary  response  is  obtained. 
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The  generalized  force  is  also  a  complex  quantity.  The  force  cross  Power  Spectral  Density 
(PSD)  is  obtained  by  integrating  the  pressure  cro^s  PSD  over  pairs  o' r..iir!  events.  ? 
gives  an  exact  formulation  of  force  cross  PSD  which  has  no  limitations  on  uppc*  response 
frequency  or  finite  element  size.  This  is  important  in  solving  the  class  of  turbulent  bound¬ 
ary  layer  problems  where  the  pressure  is  not  constant  over  the  finite  element.  01.  >n  and 
Lindbergh9)  follow  the  same  format  as  Jacobs  and  Lagerquist^)  except 4 *■> at  they  obtain 
the  force  cross  PSD  by  assuming  a  constant  pressure  over  the  finite  ele  mt  pairs.  The 
force  cross  PSD  is  obtained  by  multiplying  the  pressure  cross  PSD,  evaluated  by  using  the 
separation  distances  between  element  c.g.  ’s,  by  the  consistent  forces  for  a  unit  psi  on  the 
elements.  This  imposes  a  frequency  limit  above  which  no  valid  calculations  can  be  made. 

For  the  problem  they  analyzed,  the  upper  limit  is  1500  Hz. 

Jacobs,  et.  al.,  present  a  strong  argument  for  elimination  of  cross  modal  terms  in 

the  response  analysis.  In  Reference  16  they  show  a  variation  of  only  nine  percent  on  neak 
PSD  points  and  six  percent  on  mean  square  response  when  the  cross  modal  terms  are  ex¬ 
cluded.  In  Reference  17,  in  which  the  problem  contained  modes  very  close  in  frequency, 
they  stated  that  inclusion  of  cross  modal  terms  did  not  increase  deflection  response. 

Schvveiker  and  Davis  are  the  only  imestigalors  found  which  treat  a  class  of  problems 
where  the  pressure  PSD  (not  pressure  cross  PSD)  is  not  a  constant  over  the  entire  structure. 
They  solved  the  problem  of  a  conical  sheil  subjected  to  pressure  fluctuations  in  the  turbulent 
boundary  layer.  The  dynamic  model  is  based  on  an  assumed  closed-form  solution;  however 
the  structure  is  broken  up  into  a  series  of  sub-areas,  much  like  finite  elements,  to  formulate 
the  response.  The  normal  coordinate  of  each  sub-area  is  the  integral  of  the  assumed  mode 
shape  over  the  areas.  The  force  cross  PSD  between  pairs  of  sub  areas  is  obtained  from  the 
product  of  the  integration  of  the  spatial  correlation  function  over  the  sub-areas,  the  mto, l  o¬ 
tion  of  the  mode  shape  over  the  sub-areas  and  square  root  of  the  product  of  the  ho  :ogenocu 
power  spectra  on  the  areas,  e.g. 


. _ r  1  r.jC  <P  dA  f  v>dA^ 

/  /  c({. ,.»>«,  Mk  __ 

LaI  Ak 


where 


<l>  ,‘i'p  =  Homogeneous  power  spectra  on  sub-areas  k  and  1. 

C(|,  n,w)  =  Spatial  correlation  function. 


>p  =  Assumed  mode  shape. 


They  have  shown  the  insignificance  of  cross  modal  terms  in  the  response  formulation.  Even 
when  they  had  modes  within  1  Hz  of  each  other  a  variation  of  only  15  percent  was  observed  in 
peak  PSD  and  three  percent  in  rms  levels  by  neglecting  cross  modal  terms. 


It  becomes  obvious  from  these  observations  that  the  omission  of  iie  cross  modal  terms  from 
the  response  formulation  is  an  acceptable  simplification.  This  assumption  gi'eatly  reduces 
the  required  computation  time  and  the  computer  core  storage  thereby  resulting  in  substantial 
cost  savings  and,  as  shown  in  References  8,  15,  16  and  17,  no  appreciable  loss  in  accuracy. 
The  elimination  of  cross  modal  terms  in  the  response  technique  also  eliminates  the  imaginary 
response.  Therefore  a  co-force  cross  PSD  need  only  be  considered. 

Another  simplification  which  can  be  made  in  the  solution  of  re-entry  vehicle  response  to 
re-entry  aerodynamic  environments  is  the  assumption  of  a  constant  pressure  cross  PSD  over 
each  finite  element^) .  This  assumption  is  predicated  upon  the  high  (vs.  airplanes)  con- 

vected  velocities  (10  to  20  K  ft/sec).  The  argument  of  the  term,  cos  — — ,  in  the  spatial 

Vc 

correlation  function  does  not  change  significantly  over  a  separation  distance  equal  to  the 
length  of  the  element.  (It  should  also  be  pointed  out  that  an  upper  frequency  band  of  2000  Hz 
and  a  finite  element  length  of  approximately  three  inches  is  assumed.) 


~  =0.314  rad  or  1.8° 
vc 

The  cosine  of  18°  is  0.95,  so  no  significant  reduction  in  pressure  cross  PSD  would  occur. 


The  last  information  gained  from  the  literature  smwey  was  the  way  in  which  Reference  20 
handles  the  variable  homogeneous  pressure  PSD  over  pairs  of  finite  elements.  Their  ap¬ 
proach  is  to  use  a  pressure  PSD  between  element  pairs  equal  to  the  square  root  of  the  pro¬ 
duct  of  the  pressure  PSD  over  each  element  .  This  technique  was  implemented  in  the  re¬ 
sponse  formulations  which  follow. 


2.  OPTIMIZED  TECHNIQUE 

In  the  development  of  the  optimized  technique,  tho  simplifying  assumptions  discussed  in  the 
literature  survey  were  implemented.  These  assumptions  are: 

(1)  Cross  modal  terms  can  be  eliminated, 

(2)  Pressure  cross  PSD  is  a  constant  over  each  finite  element,  and 

(3)  The  cross  PSD  of  the  response  is  of  little  significance  and  need  not  be  computed. 

The  starting  point  in  the  response  formulation  is  the  all-inclusive  I’elationship  in  which 
there  are  no  simplifying  assumptions,  which  is 

[<I>rm]  =[*]  [n*(i«)]  [*]T  [*fm]  [*]  [h (iw) ]  (7|T  (12) 
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where 


iw 

[  H(iw)J 

f. 

M 


=  Force  cross  PSD 
=  Complex  function  of  frequency 
=  Modal  transfer  function 
=  Complex  conjugate 
-  Matrix  of  the  mode  shapes 
=  Output  cross  PSI) 


The  terms  in  the  modal  transfer  function  matrix,  [ll (i «j)]  ,  are  of  the  following  form 


(13) 


where 

w  k  =  Natural  frequency  (raa/sec) 

«  =  Calculation  frequency 

r  ••=  Modal  damping  as  a  ratio  of  the  critical  value 

M^  =  Generalized  mass,  equal  to  1  if  l  IP  1  is  normalized  with  respect  to  the 

mass  matrix . 

The  force  cross  PSD  matrix  is  obtained  from 
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where 


-  Pressure  cross  PSD,  variable  over  the  shell  element  and  a 
function  of  the  separation  distances  %  and  »/  . 


w 


n 


-  Vector  of  the  displacement  functions  used  to  obtain  the 
generalized  forces  on  the  n^1  clement. 


=  Assemblage  summation  process  used  in  setting  up  finite  element 
_ |  master  matrices  over  N  elements. 


The  force  cross  PSD  matrix  is  Henmtiun,  comprised  of  the  sum  of  a  real  or  co  matrix 
vh  ich  is  symmetric  and  an  imaginary  or  quad  matrix  which  is  skew  sy  m  mu  trie . 


(15) 


Let  us  examine  the  force  cross  PSD  in  light  of  the  simplifying  assumptions,  if  the  pressure 
cross  PSD  is  a  constant  over  euJi  finite  element  the  force  cross  PSD  can  be  expressed  as 


t!>p  (  ?»’!»  iw) 


(16) 


where 


S  >4 


=  Separation  distances  between  the  centers  of  the  two  elements 

=  Generalized  force  vectors  for  the  or  element  due  to  a 
unit  pressure  loading 


The  computer  time  savings  to  be  gained  with  this  assumption  is  that  the  integration  over 
the  pairs  of  finite  elements  car  be  eliminated.  L\ en  though  this  assumption  is  usualfy  valid 
for  re-enfcx’y  vehicles  subjected  to  aerodynamic  loading  during  re-entry ,  it  will  not  be  valid 
for  the  acoustic  chamber  test  environments.  Here  the  effective  convected  velocity  is  the 
speed  of  sound,  which,  if  a  three  inch  finite  element  length  is  assumed,  would  allow  for  an 
upper  bound  in  calculation  frequency  of  only  200  IIz.  Therefore  in  order  to  correlate  with 
test  environments,  the  form  of  force  cross  PSD  as  given  in  liquation  (14/  must  be  retained. 
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Neglecting  cross  mortal  terms  has  a  significant  effect  on  tlie  response  formulation  complexity. 
The  response  formulation  is  now  reduced  to  the  following  form: 

nmodes  T  T 

M-g  H,  tvl  I’LH,  » 

The  quad  force  cross  PSD  drops  out  because  of  its  skew  symmetry. 


Hi  [v*]  K  co 


The  above  response  formulation  greatly  reduces  the  computation  time  end  required  com¬ 
puter  core  storage  with,  as  shown  in  References  8,  15  16,  and  17,  no  significant  loss  of 
accuracy. 

The slast  simplification  results  when  Lhe  calculation  of  response  cross  PSD  is  omitted.  In 
other  words,  only  the  diagonal  elements  of  [<l>  K  (w;]  arc  computed,  onlj  in  a  few  circum¬ 
stances  would  cross  PSD  response  be  useful.  The  diagonal  terms  art,  those  used  in  test 
result  comparisons  to  assess  structural  integrity  ar.J  to  derive  component  test  specifications. 
The  response  formulation  then  becomes 


n  modes 

V">=  E 

k  =  1 


f  ^  11  2 

i  k  k 


\  M 
-  C  J 


where 


v5 ..  =  i  coordinate  in  the  k  mode  shaoe. 

ik 

‘I’s .  (W)  -  PSD  of  displacement  of  the  coordinate. 

The  different  environments  are  accommodated  through  proper  buildup  uf  the  force  cross 
PSD  matrix.  The  following  subsections  give  the  details  of  obtaining  the  force  cross  PSD 
matrix  for  the  different  environments.  Trior  to  proceeding  into  that  development,  it  would 
be  first  helpful  to  examine  the  SAJ30R  finite  element  upon  which  the  development  is  based. 

SABOR  rs  a  computer  program  developed  at  MIT  to  perform  static  anulj  set  of  shells  of  rev¬ 
olution  via  the  finite  element  method.  The  finite  element  implemented  in  the  program  is  a 
conical  frustum.  The  detailed  derivation  of  the  conical  frustum  element  is  given  in  Reference 
21.  There  are  four  generalized  coordinates  on  each  end  of  the  topical  element.  'Die  coordi¬ 
nates  and  their  sign  conventions  are  shown  in  Figure  7.  In  the  circumferential  direction  the 
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coordinates  vary  as  a  Fourier  series,  e.  g.  cos  n  0  ,  sin  n  0 ,  where  n  is  the  harmonic  wave 
number.  The  elemental  stiffness  and  mass  matrices  are  constructed  such  that  each  har¬ 
monic  is  uncoupled  from  the  other;  therefore  each  harmonic  can  be  analyzed  separately  and 
the  results  superimposed.  In  this  study  we  are  only  dealing  with  the  zeroth  harmonic,  or 
axial  model,  and  the  first  harmonic,,  oi  lateral  model.  This  is  because  the  beamlike  con¬ 
struction  of  the  component  areas  are  not  excited  by  harmonics  higher  than  one.  The 
SAI30R  shell  stiffness  is  derived  through  the  assumption  of  a  pre-described  displacement 
field.  This  same  displacement  field  is  used  in  the  derivation  of  the  generalized  forces 
under  a  distributed  loading.  This  is  accomplished  by  equating  the  virtual  work  of  the  gen¬ 
eralized  forces  at  the  element  nodes  to  the  virtual  work  of  the  distributed  loading  over  the 
element.  Reference  22  contains  the  details  of  this  formulation.  The  integrals  for  the  gen¬ 
eralized  forces  under  a  unit  pressure  loading  normal  to  the  surface  are  shown  below 
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The  remainder  of  the  forces  are  zero  for  normal  pressure  loading.  These  elemental 
forces  are  related  to  the  generalized  system  forces  by  the  transformation 


(*)*[’]  (5} 
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where 
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COS  (* 
0 

sin  « 
0 
0 
0 
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0  -sin  a  0  0 

10  0  0 

0  cos  a  0  0 

0  0  10 

000  cos  a 
0  0  0  0 

0  0  0  sin  a 

0  0  0  0 


0  0  0 

0  0  0 

0  0  0 

0  0  0 

0  -sin  a  0 
1  0  0 

0  cos  a  0 

0  0  1 


and  ot  ~  Cone  half  angle. 


(21) 


It  should  be  obsersed  that  the  circumferential  integration  is  independent  of  the  meridional 
integration.  This  is  important,  as  the  force  cross  PSD  will  be  evaluated  in  two  steps, 
first  tire  circumferential  integration  and  second  the  meridional.  Because  the  mode  shape 
variation  in  theta  is  known  (cos  n  'J  j,  this  integral  can  be  evaluated  in  closed  form  and 
multiplied  by  the  meridional  results. 

3.  RE-ENTRY  VIBRATION  RESPONSE 


In  Section  III  it  was  shown  that  the  narrowband  spatial  correlation  function  for  attached 
turbulent  flow  can  be  represented  as: 


-  ci  U1 

p  * 

C(S,  1 ,  w)  =  e  e 

In  cylindrical  coordinates  the  above  is 
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where 

5  *  and  Vc  are  now  the  average  boundary  layer  displacement  tliickness  and 
convection  velocity  between  locations  Sj  and  S2. 


(22) 
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Now  in  exact  form  the  force  cross  PSD  between  the  two  elements  k  and  m  would  be 
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(24) 
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whore 


1  /k 


-  Expressions  in  the  integrals  oi  Equations  (Tf  giver*  previous!  ■ 


v  (f)  =  Homogeneous  pressure  PSD  over  element  k, 
k 


A  typical  terra  in  the  above  relation  would  then  be 
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where  sT  is  the  meridional  distance  from  the  beginning  of  element  1  to  the  beginning 
of  element  2.  The  term  sin^  u  comes  from  the  transformation  matrix, 
Equation  <2i),  on  both  elements. 
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The  integration  with  respect  to  0  can  he  performed  separately  as 
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(26) 


Since  the  above  expression  is  independent  of  the  driving  frequency,  it  can  be  evaluated  for 
each  pair  of  elements  and  multiplied  by  the  results  of  the  meridional  integration,  which 
does  change  with  driving  frequency  .  In  performing  the  integration,  care  must  be  taken  to 
integrate  over  the  proper  areas  to  account  for  the  absolute  sign.  For  the  zeroth  harmonic 
or  n  =  0  the  integration  yields  the  circumferential  joint  acceptance: 


0 

tT  00 
j  -2 


5  *4jr 

r  c 
2  2 


"2  1 


"2ifc  r 
2  2 


-2,c2ri 


-2’°2<V1V 


+e 


-e 


■jj 

(21) 


It  is  observed  that  for  turbulent  boundary  layer  pressure  fluctuations  in  which  rA  *  is 
greater  than  one  and  c2  =  .1,  that  only  the  first  term  is  significant.  therefore  is 

approximated  as  1-2 
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Vi 


(28) 


For  higher  harmonics  in  which  n  is  1  or  greater  a  considerably  more  complex  expression 
for  ,J  ”  results. 
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2  r  c 

2  2 

'sin 

[2  7r  (n  -Hi  Y'r2)] 

sin[27r(n-n  l^/l^)]* 
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2  (n+n  r  /r2) 

2(n-nr][/r2) 

The  remainder  of  the  integration  of  Equation  (25)  is  associated  with  die  meridional  direction 
and  varies  depending  upon  the  assumption  that  the  pressure  PSD  is  constant  or  not.  If  the 
pressure  PSD  is  considered  constant  over  the  individual  elements,  the  spatial  correlation 
function  is  removed  from  the  integrals  and  evaluated  based  on  the  separation  distance 
between  the  element  centers,  e.g. 
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The  integral  expressions  can  be  evaluated  initially  for  each  of  the  generalized  forces  on  an 
element  and  storedlat  the  beginning  of  the  computations  for  use  thereafter. 

If  tire  constant  PSD  assumption  is  not  valid  (as  in  acoustic  chambers),  the  meridional 
integration  must  be  carried  out  for  each  pair  of  elements  at  each  new  driving  frequency  and 
the  spatial  correlation  function  must  of  course  be  included  in  the  integration.  The  form  of 
the  equation  thus  becomes 


The  force  cross  PSD  expressed  as  Equation  (31)  is  numerically  evaluated  using  the 
trapezoidal  rule  of  integration.  Fox  each.pair  of  elements,  sixteen  separate  integrations 
are  involved.  The  transformation  matrix  of  Equation  (21)  is  used  to  expand  the  4  x  4 
matrix  to  the  required  8x8  matrix  for  each  element  pair. 


Thus  far  only  the  force  cross  PSD  for  SABOR  finite  elements  has  been  discussed.  For  the 
beam-column  type  finite  elements  used  to  model  the  nooe  section/*  a  reduced  number  of 
coordinates  per  element  is  used.  In  the  axial  model  (zeroth  harmonic)  only  two  coordinates 
per  element  are  required.  For  the  iutcral  model  (first  harmonic)  four  coordinates  per 
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element  are  required.  If  a  generalized  force  set  were  derived  for  these  axial  and  lateral 
coordinates,  it  would  be  based  on  a  linear  displacement  function  in  the  axial  model  and  a 
cubic  displacement  function  in  the  lateral  model.  It  develops  that  the  generalized  force 
expressions  for  SABOR  radial  and  rotational  coordinates  are  identical  to  those  for  a  beam 
with  an  assumed  cubic  displacement  function.  The  SABOR  axial  generalized  force  ex¬ 
pressions  are  sufficiently  accurate  to  use  for  axial  nose  forces,  therefore  the  same 
expressions  used  for  SABOR  generalized  forces  can  be  applied  in  the  nose  sections,  elim¬ 
inating  the  need  for  separate  theory  and  checks  in  the  computer  codes. 

4.  REVERBERANT  CHAMBER  RESPONSE 


Development  for  the  force  cross-PSD  matrix  as  applied  to  an  acoustic  chamber  operating 
in  the  reverberant  mode  will  now  be  discussed.  The  spatial  correlation  function  can  be 
described  as: 


sin  —  £  sin  —  »j 

C  (?,  »,  «)  =  ~ ^ - (32) 

—  £  — v 

c  c 

where 

c  =  Sonic  velocity 
$  ,  v  =  Separation  distances 
w  =  Frequency 

In  the  cylindrical  coordinate  system  the  above  expression  would  be: 
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The  above  expression  then  replaces  the  narrowband  spatial  correlation  function  for  turbu¬ 
lent  flow  in  Equation  (25)  to  obtain  the  terms  required  in  the  force  cross  PSD  matrix  for 
reverberant  chambers.  Numerical  integration  will  be  used  meridionally,  as  in 
Equation  (31). 

Circumferentially,  the  reverberant  spatial  correlation  function  results  in  the  following  ex¬ 
pression  for  the  joint  acceptance: 
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The  closed-form  integration  of  the  above  yields  sine  and  cosine  integrals 


which  must  be  evaluated  in  series  expressions,  rational  fraction  approximations  or 
asymptotic  expansions  depending  upon  z.  Therefore  the  computer  implementation  was 
again  accomplished  via  numerical  integration. 

The  basic  response  formulation  of  Equation  (18)  then  holds  with  I"  <l>p  (u>)l  developed  as 
just  described.  L  J 

5.  PROGRESSIVE  WAVE  CHAMBER  RESPONSE 

The  narrowband  spatial  correlation  function  representative  of  an  acoustic  chamber  opera¬ 
ting  in  the  progressive  mode  can  be  described  as: 

C  (  £ ,  v  ,  «  )  =  cos  cos  (35) 


where 


i p  =  Incident  angle. 

It  is  assumed  that  the  \ehicle  is  placed  in  the  chamber  such  that  the  sound  waves  pass  over 
the  vehicle  in  tjhe  direction  of  its  axis.  Therefore  in  the  cylindrical  coordinate  system  the 
above  expression  is 
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It  is  observed  that  circumferentially  there  is  perfect  correlation.  This  means  that  the 
J  0  q  integration  of  Equation  (34)  is  simply 
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Thus  under  ideal  progressive  conditions  only  the  axial  or  zeroth  harmonic  modes  of 
vibration  can  be  excited. 
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The  force  cross  PSD  development  for  the  meridional  direction  is  the  same  as  Equation  (31), 
except  that  Equation  (3G)  is  substituted  in  place  of  the  turbulent  boundary  laye’  coirelati  ? 
functions. 

6.  COMPUTER  PROGRAMS 

The  response  formulations  developed  in  this  Section  have  been  implemented  into  computer 
programs.  A  separate  program  has  been  written  u  compute  the  re-entry  vehicle  vibration 
response  to  each  of  the  environments  discussed.  These  programs  and  associated  users 
manuals  are  available  through  the  Air  Force  Flight  Dynamics  Laboratory/F  ’ 

The  program  which  computes  the  response  of  a  re-entry  vehicle  to  re-entry  environments 
is  named  REVIB.  This  program  is  developed  to  accept  vehicle  dynamic  characteristics  of 
a  model  based  on  SABGR  shell  finite  elements  and  beam  finite  elements.  It  is  assumed  that 
the  dynamic  model  frequencies  and  mode  shapes  are  computed  using  existing  in-house 
programs  for  this  purpose  and  the  results  stored  on  tape  in  formats  compatible  with  the 
REVIB  computer  program.  In  addition  to  the  dynamic  model  input,  parameters  required  to 
describe  the  aero-acoustic  environment  must  be  input.  These  parameters ‘which  are 
described  versus  vehicle  station  include  boundary  layer  displacement  thickness,  convection 
velocity  and  overall  mean  square  dynamic  pressure.  In  addition  the  decay  constants  in  the 
spatial  correlation  functions  must  be  input.  Finally,  vehicle  and  finite  element  geometry, 
which  are  utilized  for  the  internal  generation  of  discrete  force  cross  PSD's  from  the  cli  s 
tributed  pressures,  must  be  input. 

The  input  is  required  for  both  the  forebody  and  the  aft  cover  if  both  forebody  transitional 
or  turbulent  flow  and  wake  flow  is  considered.  The  input  can  be  restricted  to  merely  the 
forebody  parameters  if  wake  flow  excitations  are  omitted  from  the  forcing  function.  In 
addition  capability  exists  to  input  discrete  force  FSDs  representative  of  shock  loadings. 

This  is  an  option  which  is  exercised  only  if  shock  loads  are  anticipated  on  the  vehicle. 

The  response  formulation  of  equation  (18)  is  computerized.  In  order  to  keep  REVIB  as 
general  as  possible  the  force  cross  PSD  formulation  described  by  equation  (14)  is  pro¬ 
grammed.  Output  from  REVIB  consists  of  response  power  spectral  density  values  at  all 
degrees  of  freedom  for  either  or  both  acceleration  and  displacement  at  all  of  the 
frequencies  selected  by  the  user  for  calculations.  In  addition,  overall  root  mean  square 
values  for  each  coordinate  power  spectral  density  are  printed.  Similar  response  output  is 
also  provided  for  the  quantities  computed  by  matrices  transfonning  acceleration  or 
displacement. 

PROG  and  REVERB  are  the  programs  which  compute  re-entry  \ehicle  vibration  response 
to  the  ground  test  acoustic  environments  associated  with  ideal  progressive  wave  chambers 
and  reverberant  chambers.  Input  to  these  programs  consists,  in  addition  to  the  finite 
element  geometry,  of  sound  pressure  levels  at  the  calculation  frequencies  for  the  vehicle 
forebody  and  its  base.  Output  is  identical  to  the  REVIB  program. 


In  addition  a  fourth  program  was  written  which  computes  the  re-entry  vehicle  vibration 
response  to  acoustic  environments  representative  of  the  Flight  Dynamics  Laboratory  Sonic 
Fatigue  Facility.  This  program  is  named  SOFAFA.  Input  is  identical  to  PROG  and 
REVERB  and  output  the  same  as  REVIB.  This  acoustic  environment  is  described  fully 
in  Section  VIII  of  the  report. 


SECTION  y 


PROTOTYPE  VEHICLE  DYNAMIC  CHARACTERISTICS 

1.  DEFINITION  OF  PROTOTYPE  VEHICLE  ' 

The  prototype  vehicle  supplied  for  this  study  was  a  Fuzing  and  Arming  Test  and  Evaluation 
Program  (FATE)  re-entry  vehicle.  It  is  a  conical  vehicle  with  a  4.5  deg”oe  half  cone  angb 
approximately  85  inches  long.  The  substructure  is  alu.  linuin  and  the  heat  shield  is  ,ai'b- 
phenolic.  There  are  three  sections  to  the  vehicle  structure  and  two  component  package- . 

The  nose  section  is  screwed  to  the  midsection  and  the  mid  ana  aft  sections  are  joined  througn 
a  bolted  flange.  One  component  package  is  essentially  shelf  traversing  almost  the  entire 
internal  midsection.  The  second  component  package  occupies  approximately  half  o.'  the  aft 
section.  More  details  and  a  schematic  of  the  vehicle  are  given  in  the  following  subsection 
(Figure  8) . 

2.  DYNAMIC  MODEL 

The  prototype  vehicle  was  modeled  using  beam  and  shell  (SABOE)  finite  elements.  Node 
identiiication  and  typical  coordinate  definition  for  axial  and  lateral  dynamic  models  are 
given  in  Figure  8.  Coordinate  identification  is  given  in  Table  I.  The  axial  dynamic  mod  Y 
has  86  degrees  of  freedom  whereas  the  lateral  dynamic  model  required  127  degrees  of 
freedom  for  its  description.  An  item-by-item  weight  breakdown  is  contained  in  Table  II. 

A  nose  fixture  was  used  for  testing  rather  than  the  prime  nose;  therefore  the  nose  fixture 
was  modeled  using  beam  finite  elements.  The  forward  frustum  is  either  completely  s  >lid 
(heavy  ballast)  or  a  very  thick  shell.  This  portion  of  the  prototype  vehicle  was  also  model-  . 
using  beam  finite  elements. 

As  seen  in  the  dynamic  model  of  Figure  8,  the  mid  and  aft  frusta  and  aft  cover  are  modeled 
using  SABOE  shell  finite  elements.  The  first  node  of  the  mid  frustum  was  transformed 
from  shell  coordinates  to  beam  coordinates  to  connect  with  the  beam-modeled  forward 
frustum.  Two  complex  joints  exist  at  the  mid- to -aft  frustum  interface  and  aft  frustum- 
to-aft  cover  interface.  Details  of  modeling  in  these  two  areas  are  shown  in  Figure  9.  Local 
discontinuties  of  these  types  usually  provide  high  flexibility,  which  considerably  affects 
gross  vehicle  bending  and  axial  modes.  The  joints  at  the  nose-to -forward  frustum  interface 
and  forward-to-mid  frustum  interface  are  good  threaded  joints  and  should  not  exhibit  high 
local  flexibility. 

The  experiment  package,  because  of  its  length  and  continuous  support  structure,  was 
modeled  using  beam  finite  elements.  The  weight  of  the  attached  components  is  distributed 
along  the  support  structure.  The  bulkhead  attaching  the  support  structure  to  the  mid-aft 
frustum  joint  (node  20)  was  modeled  using  SABOE  shell  elements,  as  shown  in  Figure  10. 

The  shell  coordinates  at  the  inner  node  are  transformed  to  beam  coordinates  to  attach  to  the 
beam  modeled  experiment  package  (node  39).  A  light  spring  was  used  in  the  lateral  model 
only,  between  the  forwa."'  end  of  the  experiment  package  and  the  forward-to-mid  fi’ustum 
joint  to  replace  the  O.OiO  inch  air  gap.  This  connection  slips  axially. 
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Figure  8.  Dynamic  Model 
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NOTE  -  v  coordinate  in  axial  model  not  used  since  this  defines  lor.-don  motion  only. 


TABLE  II.  WEIGHT  BREAKDOWN 


Item 

Weight  (Lbs . ) 

Nose,  prime 

4.07 

Nose  fixtui’e,  axial 

4.23 

Nose  fixture,  lateral 

3.69 

Forward  frustum 

45.64 

Mid  frustum 

34.50 

Aft  frustum 

37.16 

Aft  Cover 

1.65 

Experiment  Package 

36.98 

TM  Package 

24.48 

TOTAL  (PRIME)  184.48 

The  telemetry  (TM)  package  is  described  with  two  mass  points  in  the  lateral  model.  Inter¬ 
nal  components  are  mounted  on  two  connected  webs.  These  components  are  treated  as  one 
(node  41)  and  the  web  bending  stiffness  is  used  to  connect  them  to  the  package  outer  structure 
(node  40).  Components  mounted  on  the  support  bulkhead  are  also  modeled  at  node  40.  Axi¬ 
ally,  all  support  structure  is  stiff  enough  that  the  total  TM  package  is  treated  as  one  mass 
point  (node  41).  The  bulkhead  attaching  the  support  structure  to  the  aft  frustum  (node  24)  is 
modeled  using  SABOR  shell  elements,  as  shown  in  Figure  10.  The  shell  coordinates  at  the 
inner  mode  are  transformed  to  beam  coordinates  to  pick  up  the  TM  package  degrees  of  free¬ 
dom.  An  O-ring  exists  between  the  forward  end  of  the  TM  package  and  the  aft  end  of  the 
experiment  package.  This  was  modeled  as  a  linear  spring  for  the  lateral  model  and  as  a 
slip  joint  axially. 

3.  MODAL  survey  summary 

The  dynamic  model  verification  test  program  included  axial  and  lateral  modal  surveys  that 
established  vehicle  shell  and  internal  component  modal  frequencies  in  the  20  to  2000  II l  fre¬ 
quency  range.  The  mode  shapes  associated  with  resonant  frequencies  were  ascertained  in 
the  vehicle  meridional  and  circumferential  directions.  The  meridional  stations,  in  particu¬ 
lar,  were  consistent  with  the  nodes  chosen  in  the  dynamic  model  construction  to  facilitate 
ease  of  correlation.  Stations  at  either  side  of  significant  structural  or  load  transfer  discon¬ 
tinuities  were  of  particular  concern  because  experience  has  shown  that  these  areas  can  be 
large  contributors  to  poor  dynamic  model  correlation. 
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MID  FRUSTUM  TO  AFT  FRUSTUM 


•  NODES  KEPT  FOR  DYNAMIC  MODEL 
O  NODES  ELIMINATED  AFTER  JOINING  ELEMENTS 


Figure  9 .  Sabor  Modeling  of  Shell  Joints 


•  DYNAMIC  MODEL  SHELL  NODE  EXPERIMENT  PACKAGE 

®  SHELL  NODE  TRANSFORMED  TO  BEAM  COORDINATES 
O  NODES  ELIMINATED  AFTER  JOINING  ELEMENTS 


Figure  10.  Sabor  Modeling  of  Bulkhead  Attachments 
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The  modal  surveys  provided  quadrature  and  in-phase  response  versus  frequency  data  from 
which  resonant  frequencies  and  structural  damping  coefficients,  respectively,  were  deter¬ 
mined.  The  modal  surveys  also  provided  quadrature  responses  versus  vehicle  station 
(meridional  and  circumferential)  which  are  the  mode  shapes.  These  resonant  frequencies, 
mode  shapes,  and  damping  coefficients  were  used  to: 

(1)  Verify  the  dynamic  model,  thereby  providing  a  firm  basis  for  evaluating  acoustic 
effects  on  vehicle  response; 

(2)  Indicate  areas,  if  any,  where  the  modeling  techniques  can  be  improved  to  more 
closely  simulate  vehicle  characteristics;  and 

(3;  Provide  structural  damping  data  for  use  in  acoustic  response  predictions . 

The  nosetip  was  not  used  during  the  modal  surveys .  Upon  inspection  of  the  nosetip,  it  was 
concluded  that  possible  damage  or  degradation  could  occur  because  of  the  manner  in  which 
the  nosetip  is  attached  to  the  vibration  fixture.  Hence,  the  nosetip  was  replaced  by  a  GE- 
dcsigned  nosetip  adapter  with  similar  mass  and  stiffness  properties  and  allowing  simple 
attachment  to  vibration  fixture. 

The  test  was  designed  so  that  a  lightweight  shaker  armature  could  be  used.  This  weight 
was  necessarily  small  with  respect  to  the  test  specimen,  since  the  modal  frequencies  and 
mode  shapes  are  functions  of  the  entire  spring  mass  system,  which  includes  the  shaker 
armature . 

In  keeping  with  low  weight  fixtures,  the  suspension  system  for  both  axial  and  lateral  vibra¬ 
tion  consisted  of  bungee  shock  cord .  Included  in  the  lines  were  turnbuckles  for  height  ad¬ 
justment.  The  natural  frequency  of  the  suspended  system  was  less  than  5  Hz.  This  suspen¬ 
sion  system  supported  the  vehicle  in  a  vertical  position,  nose  down,  ci’eating  a  free-free 
condition  corresponding  to  free  flight.  Figure  11  shows  the  modal  survey  test  setup. 

During  the  axial  and  lateral  modal  testing,  the  location  of  one  shell  response  accelerometer 
was  changed  to  accommodate  the  various  phases  of  testing,  i.e.,  for  meridional  and  cir¬ 
cumferential  mode  shape  definition . 

Twelve  accelerometers  were  affixed  to  the  test  vehicle  for  shell  and  internal  component 
monitoring  during  both  the  modal  survey  and  acoustic  test.  Six  were  located  on  the  com¬ 
ponent  shelf  in  the  midsection  and  four  were  placed  in  the  telemetry  package  in  the  aft 
section.  Also,  one  each  was  attached  internally  to  the  midsection  and  aft  section  liners. 
Figure  12  shows  a  schematic  of  the  instrumented  vehicle. 

In  addition  to  tne  twelve  accelerometers  to  monitor  component  response,  four  strain  gages 
and  two  microphones  were  also  included  in  the  internal  instrumentation.  The  strain  gages 
were  mounted  in  die  vehicle  substructure  (Figure  12)  and  the  microphones  in  the  two  compo¬ 
nent  packages .  Neither  the  strain  gages  nor  microphones  were  utilized  during  the  modal  sur¬ 
vey;  however  they  were  used  in  the  acoustic  test  (SECTION  VII) .  Table  III  gives  the  internal 
sensor  descriptions . 
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Figure  11 .  Modal  Survey  Test  Setup 


Schematic  of  Sensor  Identification  and  Location 


TABLE  ni.  SENSOR  DESCRIPTION 


p. 


i 
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4.  FREQUENCIES,  MODE  SHAPES  AND  DAMPING 


The  modal  survey  pointed  out  the  need  for  modifying  the  dynamic  models .  The  calculated 
stiffness  coefficients  for  the  telemetry  package  bulkhead-ring  combination,  aft  field  joint 
and  aft  cover  were  modified  in  order  to  minimize  differences  between  analytical  and  ex¬ 
perimental  dynamic  characteristics  of  the  vehicle.  In  the  axial  model  the  most  sensitive 
area,  i.e.,  the  area  in  which  change  in  stiffness  had  the  most  affect  on  the  mode  shape, 
was  in  the  telemetry  package  bulkhead-ring  combination.  In  the  lateral  dynamic  model, 
the  stiffness  coefficient  representing  a  linear  spring  that  couples  the  forward  end  of  the 
experiment  package  to  the  forward  section  of  the  vehicle  was  increased  to  obtain  frequency 
similitude  between  analysis  and  test. 

Comparisons  of  the  analytically  and  experimentally  obtained  frequencies  for  the  axial  and 
lateral  models  are  presented  in  Tables  IV  and  V,  respectively,  along  with  mode  descrip¬ 
tions  and  the  experimentally  determined  structural  damping. 

Figures  13  through  1/  are  comparisons  of  the  mode  shapes  obtained  from  the  analyses  and 
the  modal  surveys.  The  modes  chosen  for  presentation  are  based  on  the  magnitude  of  the 
modal  responses  computed  b>  the  response  prediction  analyses,  i.e.,  those  modes  portray¬ 
ing  the  highest  modal  torces  and  most  significant  modal  responses  at  the  internal  acceler¬ 
ometer  locations  were  chosen. 

The  ordinate  variable  defines  relative  motion  along  the  vehicle  within  any  particular  mode. 
Thus  tne  mode  snapcs  must  be  normalized  in  some  manner.  The  normalization  technique 
used  on  the  mode  shapes,  in  the  figure  is  such  that  each  modal  mass  is  unity,  that  is 

^  M<p=  I 


where 


V  =  The  matrix  of  mode  shape  values 
M  =  The  system  mass  matrix 
I  =  The  unity  matrix 

The  following  legend  is  provided  for  interpretation  of  the  plotting  symbols  appearing  in  the 
figures . 

G  =  Quadrature  response  versus  vehicle  station  taken  during  dwelling  at 
resonant  frequency 

A  -  Quadrature  response  of  experiment  package  measurements  taken  during 
sine  sweeps 

0  >  O  -  Quadrature  response  of  telemetry  package  measurements  taken  during 
sine  sweeps 
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TABLE  IV.  COMPARISON  OF  AXIAL  EXPERIMENTAL  AND 
ANALYTICAL  FREQUENCIES 


Experimental 

(Hz) 

Analytical 

(IIz) 

Mode  Description 

Structural 

Damping 

Coefficient 

319.4 

339.0 

First  EP  Mode 

0.04 

368.0 

346.0 

First  Shell  Elastic;  TM  and  EP 

0.047 

442.0 

— 

Vei'y  Little  Motion 

0.03 

478.0 

— 

Very  Little  Motion 

0.03 

535.0 

525.0 

Aft  Cover 

0.03 

639.0 

649.0 

TM  Mode;  Shell  Elastic  and  EP 
(See  Figure  13) 

0.032 

1335.0 

1133.0 

Shell  Elastic  and  EP 

0.03 

1380.0 

1306.0 

Second  Shell  Elastic 
(See  Figure  14) 

0.037 

Note:  EP  -  Experiment  Package 
TM  -  Telemetry  Package 


The  unlabelecl  solid  line  represents  relative  displacement  associated  with  the  main  shell 
structure  of  the  vehicle.  The  solid  lines  labeled  EP  and  TM  represent  the  experiment 
package  and  telemetry  package  displacements  relative  to  themselves  and  to  the  main  struc¬ 
ture.  The  dotted  lines  are  shown  to  locate  the  vehicle  station  from  which  the  component 
packages  are  attached.  Physically,  the  length  of  the  dotted  line  represents  the  relative 
displacement  between  the  vehicle  and  the  component  shelf  attachment.  This  displacement, 
of  course,  defines  the  amount  jf  bulkhead  distortion. 

Since  the  mode  shapes  represent  the  numbers  used  in  the  response  analysis  to  describe 
coordinate  displacements,  etc,,  those  modes  which  exhibit  the  largest  shell  deviations  from 
the  ox’igin  of  the  ordinate  axis  would  be  expected  to  show  the  greatest  response.  This  is 
especial''  true  lAr  the  last  30  inches  of  the  vehicle  as  the  greatest  surface  area  is  present, 
thus  res.  ng  in  larger  forces. 
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TABLE  V.  COMPARISON  OF  LATERAL  EXPERIMENTAL  AND 
ANALYTICAL  FREQUENCIES 


Experimental 

m 

Analytical 

m 

Mode  Description 

Structural 

Damping 

Coefficient 

74 

75 

EP  Rocking 

0.068 

— 

98 

TM  Motion 

— 

135 

134 

Vehicle  First  Elastic 
(See  Figure  15) 

0.059 

to 

►£- 

O 

232 

Vehicle  Second  Elastic 

0.029 

260 

268 

EP  Motion;  Vehicle  Elastic 

0.03 

326 

353 

EP  First  Elastic 
(See  Figure  16) 

0.045 

460 

TM  and  Vehicle  Elastic 

0.04 

534 

580 

Vehicle  Third  Elastic 
(See  Figure  17) 

0.026 

730 

731 

EP  Second  Elastic;  TM  and 
Vehicle  Elastic 

0.0  74 

926 

1034 

Vehicle  Fourth  Elastic 

0.055 

— 

1118 

Aft  Cover  (Not  Surveyed) 

-- 

1207 

1228 

Vehicle  Fifth  Elastic,  EP 

0.038 

1425 

1447 

Vehicle  Elastic,  EP 

0.053 

Note:  EP  -  Experiment  Package 
TM  -  Telemetry  Package 
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VEHICLE  STATION.  INCHES 


Figure  13 .  Axial  Mode  Shape 


VEHICLE  STATION,  INCHES 


Figure  14.  Axial  Mode  Shape 


-4 


Figure  17.  Lateral  Mode  Shape 

In  Table  IV,  it  is  seen  that  axial  experimental  frequencies  442  Hz  and  478  Hz  have  no  cor¬ 
responding  analytical  frequencies,  and  that  their  mode  shapes  portray  very  little  motion.  It 
was  not  fully  understood  why  the  above  frequencies  were  not  predicted,  since  there  was  ex¬ 
cellent  frequency  similitude  in  the  other  six  modes.  However,  contribution  from  the  442  and 
478  Ilz  modes  will  be  negligible.  This  was  verified  in  the  acoustic  test.  That  is,  since 
there  is  very  little  motion  in  the  modes,  modal  forces  and  modal  displacements  will  be 
diminutive . 

In  Table  V,  there  is  no  corresponding  lateral  experimental  frequency  for  the  98  Hz  analytical 
prediction.  The  analytical  mode  shape  showed  very  little  motion  of  the  vehicle  nose.  Ilcnce, 
the  shaker  located  at  the  vehicle  nose  during  modal  survey  testing  would  not  produce  large 
telemetry  package  (TM)  motion.  The  only  way  to  duplicate  this  mode  experimentally  would 
be  to  physically  attach  the  shaker  to  the  TM.  Also,  in  Table  V,  the  experimental  frequency 
of  460  Hz  does  not  have  a  corresponding  analytical  frequency. 

The  structural  damping  v  alues  shown  in  Tables  IV  and  V  are  obtained  from  the  in-phase  or 
co  response  peaks  just  before  and  just  after  resonance.  The  co  response  is  zero  at  resonace. 


Structural  damping,  g',  is  related  to  peak  co  response  frequencies,  wa  and  by 


Volume  II  of  this  report  contains  the  in-phase  response  plots  from  which  the  structural 
damping  coefficients  in  the  tables  were  obtained. 
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SECTION  VI 


RE-ENTRY  VIBRATION  RESPONSE  PREDICTIONS 

1.  COMBINED  FOREBODY  AND  WAKE  ENVIRONMENTS 

The  prototype  vehicle  dynamic  model  was  subjected  to  the  combined  turbulent  boundary  lajer 
and  wake  environments  ol  Section  III.  The  response  formulation  of  Section  IV  utilizing  the 
constant  pressure  cross  PSD  assumption  was  followed.  Figures  18  through  30  show  the 
resulting  predicted  acceleration  power  spectral  densities  for  each  internal  accelerometer 
location.  Shown  on  the  figures  in  addition  to  the  PSD's  are  the  rms  accelerations.  The 
predicted  strain  PSD's  were  not  plotted,  because  they  were  too  small,  less  than  lO"1*’ 

(in/in)  2/lIz. 

In  order  to  discuss  the  predicted  responses,  it  would  first  be  helpful  to  present  the  trans¬ 
formations  utilized  in  obtaining  them.  The  dynamic  model  does  not,  in  all  instances,  have 
coordinates  which  directly  correspond  to  accelerometer  locations.  For  this  reason, 
geometric  transformations  have  been  used  to  predict  the  response  at  the  exact  accelero¬ 
meter  locations.  Table  VI presents  the  transformations  used,  where  the  qj  coordinates 
represent  degrees  of  freedom  in  the  dynamic  models . 


10  100  1000  HUS  0  42  G'S  10000 

FREQUENCY  III 


Figure  18.  Accelerometer  1AL  Response 
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ACCELERATION  SPECTRAL  DENSITY.  G"2/Hz 


FREQUENCY.  Hj 


Figure  19.  Accelerometer  1AT  Response 


.  Accelerometer  3AL  Response 
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Figure  20 


ACCELERATION  SPECTRAL  DENSITY.  G--2/HJ! 


Figure  22.  Accelerometer  4AL  Response 


ACCELERATION  SPECTRAL  DENSITY.  G”2/Hz 


ACCELERATION  SPECTRAL  DENSITY  G”2/Hz  ACCELERATION  SPECTRAL  DENSITY.  G**2/Hz 


TABLE  VI.  TRANSFORMATIONS  FOR  ACCELEROMETER 
RESPONSE  PREDICTIONS 


Accelerometer 

Axial  Model 
Coordinates 

Lateral  Model 

Coordinates 

1AL 

q79 

2-°“ll2 

1AT 

- 

0.47  qm+ 0.53 

3AL 

q83 

°-625  qi20 

3  AT 

- 

0.38  qm  + 0.62  qu9 

4AL 

q85 

0.625  q^ 

4  AT 

•• 

0.40  W  0.51  qi23 

5AL 

q36 

- 

5  AT 

- 

qi27 

GAL 

q86 

3-°  “126 

6  AT 

- 

“125  -  6 -128  “126 

8AN 

Q35 

q53 

9AN 

q56 

q8l 

10  AL 

q76 

qi07 

Table  VI  shows,  for  example,  that  accelerometer  1AL  represents  the  sum  of  the  accelera¬ 
tions  at  coordinate  79  in  the  axial  model  and  a  2  inch  radius  times  the  rotation  of  coordinate 
112  in  the  lateral  model.  Accelerometer  1AT  receives  no  contributions  from  the  axial  model 
but  is  physically  located  between  coordinates  111  and  113,  so  that  the  transformation  repre¬ 
sents  an  averaging  of  the  accelerations  at  these  coordinates.  The  rest  of  the  accelerometer 
locations  are  similarly  obtained.  It  is  observed  that  accelerometers  5AL,  and  6AL  are 
functions  of  axial  model  accelerations  only  and  1AT,  3 AT,  4AT,  5 AT  and  SAT  are  functions 
of  the  lateral  model  accelerations  only.  These  results  are  reflected  in  the  response  pre¬ 
dictions  shown  in  Figures  19  through  30,  where  it  is  seen  that  PSD  peaks  of  Figures  24  and 
26  are  due  to  only  axial  model  natural  frequencies,  whereas  Figures  19,  21,  23,  25,  and  27 
reflect  lateral  model  frequencies  at  the  peaks  of  the  PSD's.  Figures  18,  20,  21,  28,  and  29 
contain  peaks  associated  with  both  axial  and  lateral  natural  frequencies,  since  these  loca¬ 
tions  are  coupled  in  the  axial  and  lateral  directions.  Figure  30,  which  represents  the 
accelerometer  located  on  the  aft  cover,  only  exhibits  one  dominant  frequency,  which  is  the 


aft  cover  frequency.  This  is  expected,  since  the  aft  cover  mode  at  524  Hz  in  the  axial  model 
responds  so  significantly  that  it  masks  the  contributions  of  other  modes  at  their  natural  fre¬ 
quencies.  If  there  were  no  aft  cover  loading,  Figure  30  would  exhibit  peaks  at  other  modes 
in  addition  to  the  one  at  524  Hz. 

2.  SEPARATE  FOREBODY,  WAKE  ENVIRONMENTS 

The  responses  shown  in  the  figures  are  for  a  simultaneous  application  of  the  re-entry  en¬ 
vironments  in  Section  III;  that  is,  the  turbulent  boundary  layer  environment  was  applied  at 
the  same  time  as  the  wake  environment.  The  effects  of  each  environment  independently  ap¬ 
plied  will  now  be  discussed. 

Basically  little  difference  would  be  expected  in  the  results  for  Figures  18  through  29  if  only 
the  forebody  environment  were  applied.  The  aft  cover  environment,  as  shown  in  Section  III, 
is  significantly  lower  than  the  forebody  environment.  The  force  cross  PSD  matrix  for  the 
combined  environments  would  have  smaller  terms  for  the  part  of  the  matrix  associated  with 
the  aft  cover  coordinates  versus  the  terms  associated  with  the  forebody.  Also,  except  for 
analytical  modes  identified  as  aft  cover  modes,  littie  aft  cover  motion  exists  in  the  other 
predicted  modes .  Thus  when  the  modal  force  is  obtained  through  the  pre-  and  post¬ 
multiplication  of  the  force  cross  PSD  matrix  by  the  individual  mode  shapes,  there  would  be 
little  contribution  from  the  wake  environment.  The  responses  predicted  for  the  combined 
environments  are  therefore  essentially  the  same  as  would  be  expected  for  the  attached  tur¬ 
bulent  flow  alone  for  the  internal  component  accelerometers. 

This  is  obviously  not  the  case  for  the  aft  cover  accelerometer.  Under  the  turbulent  bound¬ 
ary  layer  environment  alone,  the  response  predictions  for  the  aft  cover  would  not  show  the 
dominance  of  the  aft  cover  mode .  The  response  levels  would  be  greatly  reduced  over  that 
shown  in  Figure  30,  and  the  PSD  would  be  characterized  by  peaks  at  each  natural  frequency 
whose  mode  shape  has  oven  slight  aft  cover  motion. 

Under  the  reverse  situation,  where  the  environment  is  only  applied  to  tie  base  of  the  re¬ 
entry  vehicle,  the  opposite  trends  would  result.  Figure  30  would  be  unchanged,  but  Fig¬ 
ures  18  through  29  would  be  reduced  by  orders  of  magnitude  and  reflect  modes  which  contain 
a  slight  amount  of  aft  cover  motion . 
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SECTION  VII 


RESPONSE  VERIFICATION  TEST 


1.  OBJECTIVES 

The  response  verification  test  was  an  acoustic  test  utilizing  the  AFFDL  Sonic  Fatigue 
Facility. 

The  objectives  were  to: 

(1)  Verify  the  optimized  response  prediction  technique  by  obtaining  vehicle  response 
data  plus  internal  vibration  response  in  a  form  suitable  for  correlation  with 
analyses;  and  to 

(2)  Define  vehicle  simulation  techniques  by  making  accurate  measurements  of  external 
and  internal  acoustic  levels  and  spectra  in  order  to  define  and  express  in  sufficient 
detail  the  experimental  forcing  function. 

2.  TEST  REQUIREMENTS 

a.  Acoustic  Test  -  AFFDL  Sonic  Fatigue  Facility 

The  vehicle  was  subjected  to  two  separate  test  environments  of  reverberant  and 
progressive  wave  sound  fields  in  the  Large  Test  Chamber  at  the  AFFDL  Sonic  Fatigue 
Facility.  The  support  system  consisted  of  three  lengths  of  bungee  cord  and  three  turn- 
buckles.  One  bungee  cord  was  attached  to  the  nose  tip  and  the  two  remaining  to  the  aft  end. 
Figure  31  shows  the  support  system.  The  vehicle  was  positioned  in  the  chamber  with  its 
nose  approximately  10  feet  from  the  siren  bank.  Its  roll  axis  was  coincident  with  the  geo¬ 
metric  center  of  the  sirens. 

hi  addition  to  the  internally  installed  instrumentation  to  monitor  vehicle  response, 
airays  of  microphones  were  utilized  at  the  external  vehicle  surface  to  obtain  the  acoustic 
environment  description.  Figure  32  shows  the  external  instrumentation.  Based  on  sound 
pressure  data  from  the  test  facility  and  response  predictions  for  the  acoustic  environments, 
a  minimum  sound  pressure  level  spectrum  was  defined  for  the  verification  test.  This  is 
shown  in  Figure  33.  It  was  also  shown  in  the  pretest  predictions  that  if  the  same  sound 
pressure  levels  were  applied  to  the  base  of  the  re-entry  vehicle  as  the  forebody,  the  aft 
cover  response  would  bo  in  excess  of  the  re-entry  vibration  predictions  for  wake.  A  baffle 
was  therefore  constructed  to  reduce  the  incident  sound  pressure  levels  on  the  aft  cover  in 
the  acoustic  test  by  15  to  18  dB  over  that  on  the  forebody. 

b.  Chamber  Environments  -  Reverberant,  Progressive  Wave 

In  order  to  achieve  a  diffuse  sound  field,  the  AFFDL  Large  Test  Chamber  curtains 
were  raised  and  the  floor  uncovered.  In  the  progressive  test  the  curtains  were  down  and  the 
floor  was  covered.  The  same  sound  piessure  spectra  were  required  for  each  test;  they  are 
given  in  Figure  33. 
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Figure  32.  Schematic  of  External  Acoustic  S 
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Figure  33.  Acoustic  Spectrum  —  Minimum  Level 


The  data  processing  requirements  for  each  test  environment  were  as  follows:  All 
raw  acoustic  analog  data  were  analyzed  to  obtain  Sound  Pressure  Level  (SPL)  on  a  1/3-ootavc 
band  basis  for  post  test  quick-look  evaluation.  Power  spectral  density  (PSD)  analyses  were 
later  performed  on  all  acoustic  sensor  data.  Figure  34  shows  quick-look  acoustic  data 
processing  requirements.  In  addition,  detailed  data  processing  of  the  acoustic  sensors 
comprising  the  L-shaped  patch  arrays  were  performed.  The  intent  was  to  determine  spatial 
correlation  functions  and  cross- PSDs  in  order  to  express  the  experimental  forcing  function 
in  a  form  suitable  for  use  with  the  FATE  HI  A  dynamic  model.  This  was  the  primary 
reason  for  recording  the  L-shaped  array  microphones  on  specific  tape  recorder  channels. 
Figure  35  shows  detailed  acoustic  data  processing  requirements. 


Figure  34.  Quick- Look  Acoustic  Data  Processing 


NOTE:  NARROW  BAND  SPATIAL  CORRELATION  COEFFICIENTS  ARE  TO  BE 
OBTAINED  AT  THE  FOLLOWING  FREQUENCIES: 

125,  500,  1000  AND  2000  Hz 


Figure  35.  Detailed  Acoustic  Data  Processing 

Since  acoustic  sensors  LI  and  L8  were  the  pivot  sensors  in  the  L-shaped  patch 
arrays,  the  following  combinations  were  analyzed: 

NO.  1  PATCH  ARRAY 


Meridional 


Circumferential 


LI  -  L2 
LI  -  L3 
LI  -  L4 
LI  -  L5 


LI  -  L6 
LI  -  L7 


NO,  2  PATCH  ARRAY 


Meridional 


Circumferential 


L8  -  L9 
L8  -  L10 
L8  -  Lll 
L8  -  L12 


L8  -  L13 
L8  -  L14 
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c.  Vibration  Response 


The  prototype  vehicle  response  to  the  acoustic  environments  was  monitored  via  the 
internally  mounted  accelerometers  and  strain  gages.  The  acceleration  and  strain  data,  alter 
initial  RMS  determination,  were  analyzed  for  frequency  content  by  obtaining  power  spectral 
densities.  Figure  36  shows  the  data  processing  requirements  for  the  acceleration  and  strain 
data. 

3.  TEST  RESULTS 

a.  Acoustic  Environment 

Upon  completion  of  the  data  evaluation  from  the  two  acoustic  tests,  it  became  clear 
that  there  was  little  difference  in  the  environment  between  the  two  tests.  The  specimen  was 
mounted  too  close  to  the  siren  bank  to  achieve  a  truly  diffuse  sound  field.  The  two  test 
environments  can  be  best  described  as  modified  progressive.  The  detailed  acoustic  data 
evaluation  is  given  in  Section  VIII.  Volume  II  of  this  report  contains  the  entire  data  obtained 
during  the  acoustic  testing. 

Figures  37  and  38  represent  typical  one-third  octave  band  sound  pressure  level 
spectra  obtained  during  each  test  from  a  forebody  microphone.  Figures  39  and  40  represent 
sound  pressure  spectra  obtained  on  the  aft  cover  during  each  test.  It  is  seen  that  the  baffle 
performed  as  designed. 

Attention  will  now  be  turned  to  the  internal  acoustic  environments.  The  purpose  of 
microphones  2PL  and  7PL,  mounted  in  the  component  packages,  was  to  ascertain  the  amount 
of  the  acoustic  environment  that  passes  through  the  shell  structure.  This  environment  could 


NOTES:  (1)  DATA  ANALYZED  FROM  20  Hz  to  2  kHz 
(2)  FILTER  REQUIREMENTS 

20  TO  200  Hz  WITH  5  Hz  FILTER 
200  Hz  TO  2  kHz  WITH  20  Hz  FILTER 

Figure  36.  Vibration  Data  Processing 
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Figure  37.  Forebody  Sound  Pressure  Leva!  In  Progressive  Test 


Figure  :J8.  Forebody  Sound  Pressure  Level  Li  Reverbei’ant  Test 


Figure  39.  Base  Sound  Pressure  Level  In  Progressive  Test 


t«J0 


Figure  40.  Base  Sound  Pressure  Level  In  Reverberant  Test 


be  of  high  enough  magnitude  to  directly  excite  the  component  packages.  Figures  41  and  42 
show  octavo  band  sound  pressure  levels  for  the  two  Internal  "librophones 

during  the  progressive  acoustic  test.  It  is  seen  that  there  is  significant  reduction  in  the 
IZs  S  external  environment  of  Fibres  37  and  38  The  Overall  Sound  Pressure 
Level  is  reduced  by  36  dB,  with  a  corresponding  29  to  32  dB  reduction  in  the  one-thiid 

octave  bands. 

b.  Vibration  Hesoonso 

As  would  be  expected  when  the  environments  are  essentially  the  same,  the  vehicle 
response  was  nearly  identical  for  both  tests.  Acceleration  power  spectra  densities  Vt ex e 
obtained  for  each  accelerometer  location.  The  figures  in  Section  VIH  contain  these  1 SD 
As  anticipated,  the  strain  gage  response  was  negligible,  so  no  attempt  was  made  to 

this  data. 


Figure  41.  Internal  Sound  Pressure  Level  (2 PL)  In  Progressive  Test 
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SECTION  VIII 


PREDICTION  TECHNIQUE  VERIFICATION 


1.  TEST  ACOUSTIC  ENVIRONMENT  DESCRIPTION 


The  acoustic  test  environment  has  been  defined  in  a  form  suitable  for  response  predictions. 
The  two  bits  of  data  required  to  make  response  predictions  are  the  pressure  power  spectral 
densities  and  the  narrow  band  spatial  correlation  functions.  The  pressure  PSDs  are  of 
course  available  di recti}  from  the  processed  data.  Figure  -13  shows  a  typical  pressure  PSD 
obtained  using  a  narrow  band  filter  for  both  the  progressive  and  reverberant  acoustic  tests. 
It  is  observed  that  neither  PSD  represents  a  smooth  curve.  Each  of  the  peaks  and  valleys 
must  be  accounted  for  in  order  to  obtain  accurate  correlations.  It  is  also  seen  that  little 
difference  exists  between  the  progressive  and  reverberant  tests . 


Attention  will  now  be  turned  to  the  derivation  of  the  narrow  band  space  correlation  functions. 
The  narrow  band  space-time  correlation  data  was  used  for  this  purpose.  Under  ideal  pro¬ 
gressive  wave  conditions  the  narrow  band  space -time  correlation  functions  should  be  of  the 
following  form  in  the  meridional  direction: 


i<  ftiu.r) 


sin  [  2  ff  (r  -S/g)] 
2nff  (r-$/2) 


cos 


2  n  fc  (  t  -  $  /  c  ) 


(39) 


where 

fj-  =  Filter  bandwidth  used  in  the  data  processing 
fc  =  Center  frequency 

$  =  Separation  distance  between  the  two  microphones 
~c  =  Sonic  velocity. 

The  above  expression  peaks  at  1.0  when  T  =  S/c'.  Figure  44  shows  the  test  space-time 
correlation  functions  at  1000  Hz  for  microphones  Ll  and  L3,  and  for  LI  and  L5.  Also  shown 
are  points  representing  the  above  expression.  The  filter  bandwidth  used  in  the  evaluation 
was  l/3  octave.  It  is  seen  on  Figure  44  that  the  test  points  very  closely  coincide  with  the 
theoretically  perfect  progressive  points.  The  only  difference  is  that  the  test  data  does  not 
peak  at  1.0  at  values  of  T  -  % /c  ,  but  is  somewhat  decayed.  This  suggests  a  narrow  band 
space  correlation  function  of  the  ideal  progressive  form,  but  with  a  damping  term  ipclnrW) 
This  damping  term  can  be  obtained  from  the  values  of  the  space-time  correlation  data  at 
r  =  |  /c  .  A  plot  of  the  values  of  the  space-time  correlation  data  at  t  -  £  /c  versus  £u>/e 


Cl£<o) 


is  constructed.  A  curve  is  then  fit  through  these  points.  This  yields  a  narrow  band  space 
correlation  function  offjjrie  following  form  in  the  meridional  direction: 

-.0364|  ?! <«> 

C($,«)  =  e  C  cos  (40) 

A  check  on  the  above  expression  is  obtained  from  the  space-time  correlation  data  at  r  -  o. 
The  narrow  band  space  correlation  function  is  the  space  time  correlation  at  *  -  0 .  Figure 
45  shows  the  final  meridional  space  correlation  function  with  the  test  data  at  r  -  o  super¬ 
imposed.  Also  shown  are  the  points  used  to  arrive  at  the  damping  term.  It  is  seen  that 
nearly  all  of  the  test  points  fall  on  the  analytical  curve. 

Under  perfect  progressive  conditions  the  circumferential  space-time  correlation  data  should 
be  of  the  following  form: 


sin  [2/rf^  t] 


2  4  6  8  in  1?  14  ;b 

iu 


C 

Figure  45.  Meridional  Spatial  Coi’relation  Function 
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Notice  the  absence  of  any  dependence  on  the  microphone  separation  distance .  The  above 
expression  is  simply  the  auto  correlation  function.  It  is  seen  that  this  expression  also 
peaks  at  1.0,  but  at  r  =  0.  Figure  46  shows  circumferential  space-time  correlation  data 
at  1000  Hz.  It  is  seen  that  the  peak  does  indeed  occur  at  t  =  0  but  it  is  not  1.0 .  This  again 
suggests  damped  circumferential  space  correlation  function.  Figure  47  shows  the  decay 
envelope  for  the  circumferential  direction  and  the  test  points  used  to  obtain  it.  The  analyt¬ 
ical  expression  is 


-0.0705  gM_ 

C(>J,w)  =  e  ’  (42) 

The  total  narrow  band  spatial  correlation  function  for  the  progressive  acoustic  test  there¬ 
fore  is 


-0.0364  <oUl  -0.0705 

C(S,n,w)=e  C  o  0  cos  . £  ^  (43) 

c 

It  is  this  expression  which  will  be  used  in  the  correlation  analysis. 

The  reverberant  test  data  is  very  similar  to  the  progressive  data.  As  shown  in  Figure  43 
the  pressure  PSD's  are  nearly  the  same.  Figures  48  and  49  show  a  comparison  between 
meridional  and  circumferential  space-time  correlation  functions  for  the  progressive  and 
reverberant  tests.  The  only  significant  difference  between  the  two  sets  of  test  data  is  the 
larger  amount  of  decay  ir.  the  peaks  that  the  reverberant  test  exhibits .  This  suggests  that 
the  reverberant  test  was  in  reality  a  more  decayed  progressive  test.  This  result  is  ex¬ 
pected,  due  to  the  close  proximity  of  the  re-entry  vehicle  to  the  siren  bank.  Another  fact 
which  verifies  that  the  reverberant  test  is  actually  a  progressive  test  is  the  peak  of  the  cir¬ 
cumferential  space-time  correlation  data  at  r  -  0 .  If  the  environment  were  reverberant, 
peaks  at  values  of  r  greater  than  zero  would  be  expected.  This  is  because  in  a  reverberant 
environment  the  microphone  separation  distance  influences  the  space-time  correlation  peaks 
in  all  directions,  not  just  the  meridional  as  in  the  progressive  test. 

The  reverberant  space-time  correlation  data  was  similarly  analyzed  as  the  progressive  to 
yield  a  narrow  band  spatial  correlation  function  of  the  following  form: 

0495 lily  -0.126  M.y. 
c  c 

C  (%,  i),  w)  =  e  e  cos  i£_L 

Thc  only  difference  between  the  above  expression  and  the  progressive  expression  is  the 
decay  constants. 


Meridional  Space- 


2.  TEST  VIBRATION  RESPONSE  PREDICTION 


Several  modifications  had  to  be  made  to  the  existing  response  prediction  techniques  in  order 
to  have  the  capability  to  predict  the  test  response.  These  modifications  were  made  .to  the 
existing  progressive  wave  response  formulation  and  were  associated  with  the  incorporation 
of  the  acoustic  test  spatial  correlation  functions. 

The  difference  between  the  acoustic  test  environment  and  the  perfect  progressive  wave 
environment  is  the  exponential  decay  term  in  both  the  meridional  direction  and  circumferen¬ 
tial  direction.  The  prediction  technique  revision  in  the  meridional  direction  was  easily  im¬ 
plemented,  as  it  just  entailed  the  addition  of  this  term  to  the  correlation  function.  In  the 
circumferential  direction,  the  change  was  more  significant.  Fortunately  the  turbulent 
boundary  layer  (TBL)  environment  also  contains  an  exponential  decay  term  so  the  closed 
form  solutions  of  the  circumferential  contributions  to  the  force  cross  PSD  for  turbulent 
boundary  layer  are  applicable  for  the  acoustic  test  with  small  alterations*  The  turbulent 
boundary  layer  expression  for  the  circumferential  spatial  correlation  function  is 


C(j?,w)  =  e  5  * 
whereas  for  the  acoustic  test 

C(v,u)  =e~C21J?4^ 

It  is  seen  that  if  C£ /(,  *  is  replaced  by  cgKo/c,  the  closed  form  expressions  for  turbulent 
boundary  layer  presented  in  Section  IV  are  applicable.  These  expressions  have  been 
incorporated  into  the  computer-implemented  progressive  response  formulation  to  result  in 
a  program  with  the  capability  of  predicting  the  acoustic  test  response. 

This  program,  coupled  with  an  average  of  the  pressure  PSD  of  all  the  microphones  irom  the 
progressive  test,  was  used  to  predict  the  test  response. 

3.  TEST-ANALYSIS  CORRELATION 

The  acoustic  test  predictions  were  done  only  for  the  progressive  test.  Generally,  good 
correlation  was  obtained  between  test  data  and  the  analytical  predictions.  There  are  a  few 
instances  where  the  correlation  was  poor,  but  these  are  readily  explainable.  Figures  50 
throughout  62  show  the  test  results  and  analytical  predictions .  The  analytical  results  were 
obtained  using  the  techniques  of  the  previous  subsection.  In  some  cases  the  analytical  mode 
shapes  were  altered  to  coincide  with  the  test  modes  from  the  modal  survey .  Each  of  the 
figures  will  now  be  discussed. 
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Figure  54 .  4AL  Vibration  Response, 
Test  Vs.  Analysis 


Figure  56.  5A'L  Vibration  Response, 
Test  Vs.  Analysis 
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Figure  55.  4AT  Vibration  Response; 
Test  Vs.  Analysis 


Figure  57.  5 AT  Vibration  Response, 
Test  Vs.  Analysis 
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Figure  58.  6AL  Vibration  Response, 
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Figure  59.  6 AT  Vibration  Response 
Test  Vs .  Analysis 
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Figure  60.  8AN  Vibration  Response, 
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Figure  61.  9 AN  Vibration  Response 
Test  Vs.  Analysis 
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Figure  62.  10AL  Vibration  Response,  Test  Vs.  Analysis 

Figures  50  and  51  show  analytical  predictions  for  accelerometers  1AL  and  1AT  along  with 
the  test  data.  It  is  seen  that  good  correlation  between  test  peaks  and  analytical  peaks  exists 
below  1000  Hz.  There  are  some  test  peaks  which  do  not  have  a  corresponding  analytical 
prediction.  One  such  frequency  is  460  IIz.  This  test  mode  was  obtained  in  the  modal  survey 
but  no  corresponding  mode  was  obtained  in  the  dynamic  model.  At  the  time  of  dynamic 
model  correlation,  it  was  believed  that  the  460  Hz  mode  would  not  exhibit  large  responses 
in  the  acoustic  test.  The  test  has  proved  this  to  be  incorrect.  On  1AL  a  test  peak  at  850  Hz 
is  noted.  This  frequency  did  not  show  up  in  the  modal  survey,  and  is  likely  a  frequency 
associated  with  the  other  orthogonal  lateral  direction.  Between  1000  and  2000  IIz  more 
discrepancy  is  noted  between  predictions  and  test  data.  This  is  to  be  expected,  because  the 
analytical  modes  were  not  as  well  correlated  with  modal  survey  test  modes.  In  fact,  above 
1500  IIz  no  test  modes  were  found  in  the  modal  survey,  but  there  are  modes  evident  in  the 
acoustic  test.  These  modes  could  possibly  be  higher  order  component  modes  which  were 
not  excited  by  the  re-entry  vehicle  nose  input  used  in  the  modal  survey.  The  trend  of  analyt¬ 
ical  predictions  being  less  than  test  data  above  1500  Hz  is  present  on  almost  all  figures,  and 
is  attributed  to  these  modes. 

Figures  52  and  53  show  analytical  and  test  data  for  accelc  cemeters  3AL  and  3 AT.  Again 
good  correlation  is  noted  at  frequencies  below  1500  Hz.  On  these  two  figures  the  460  IIz 
mode  is  more  evident.  The  dropoff  above  1500  Hz  is  also  noted. 
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Figures  54  and  55  show  comparisons  for  accelerometers  4AL  and  4AT.  Poor  correlation 
is  noted  for  4AL,  It  is  also  observed  that  the  test  peaks  for  4AL  occur  primarily  at  lateral 
frequencies:  268,  579,  750,  950,  and  1050  IIz.  The  lateral  contribution  in  this  axial  ac¬ 
celerometer  is  realized  analytically  by  multiplying  the  rotational  coordinate  in  the  lateral 
model  at  this  accelerometer  location  by  a  moment  arm  (0.625  inches)  to  account  for  the  off¬ 
set  from  the  EP  axis  (see  Section  VI) .  In  order  to  significantly  raise  the  response  at  these 
lateral  frequencies,  the  analytical  i’otational  coordinate  would  have  to  be  increased  in  each 
mode  or  the  value  of  4AL  due  to  lateral  motion  known.  There  is  no  way  to  know  exactly  how 
much  to  alter  this  rotational  coordinate  based  on  the  modal  survey,  since  it  is  difficult  to 
measure  fhis  quantity,  and  only  lateral  motion  at  the  three  ISP  stations  were  defined  in  the 
lateral  modal  survey.  This  was  one  of  the  shortcomings  of  the  modal  survey,  and  will  be 
discussed  more  fully  in  Section  XI. 

Figures  56,  57,  58,  and  59  show  the  correlation  for  accelerometers  5AL,  5 AT,  6AL,  and 
6 AT,  respectively.  These  accelerometers  are  located  in  the  TM  package.  It  is  on  these 
plots  that,  in  some  cases,  the  test  mode  shape  values  were  substituted  into  the  analytical 
mode  shape  description  to  give  better  correlation  with  the  modal  survey  data.  It  is  seen 
that  good  correlation  exists  on  Figures  58  and  59  and  poor  correlation  on  Figure  56  and  57. 
On  all  of  the  figures  the  peak  at  460  Hz  is  noted,  which  again  is  not  analytically  predicted. 

It  is  also  noted  that  the  axial  accelerometers  exhibit  peaks  at  lateral  frequencies.  There 
is  a  rotational  coordinate  in  the  dynamic  model  which  corresponds  to  the  6AL  and  6AT  ac¬ 
celerometer  locations.  This  coordinate  was  therefore  used  in  conjunction  with  the  proper 
moment  arm  to  include  in  the  6AL  response  the  lateral  frequencies.  Unfortunately,  no  such 
coordinate  in  the  dynamic  modal  was  included  for  the  SAL  and  SAT  accelerometers;  there¬ 
fore  5AL  does  not  analytically  exhibit  lateral  frequency  response  and  so  does  not  correlate 
well  with  the  test  data. 

Figure  63  shows  a  comparison  between  the  analytical  mode  shape  and  experimental  one  at 
650  Hz.  It  is  observed  that  the  analytical  mode  shape  predicts  the  same  response  for  5AJL 
and  6AL  which  is  between  the  two  test  values.  These  two  test  values  were  substituted  within 
the  analytical  mode  to  improve  the  correlation.  This  was  one  example  of  substituting  test 
values.  It  was  in  the  TM  area  where  the  substitution  of  the  test  values  was  requix’ed.  This 
was  due  to  the  inherent  complexity  associated  with  modeling  the  TM  package .  An  interesting 
observation  can  be  made  by  comparing  5AT  test  response  at  730  Hz  with  the  test  mode  shape 
at  that  frequency  in  Volume  II  of  this  report.  It  is  seen  that  5AT  exhibits  no  response  in  the 
modal  survey  at  730  Hz,  yet  responds  significantly  in  the  acoustic  test.  No  doubt  the  dif¬ 
ference  in  the  two  environments  accounts  for  this  phenomena.  It  is  also  seen  that  neither 
6AL  or  6AT  exhibit  any  response  at  all  at  2000  Hz  in  the  modal  survey,  yet  have  large  peaks 
in  the  acoustic  test.  This  pi’oblem  in  achieving  good  correlation  in  the  TM  package  points 
out  the  need  for  modal  suxweys  of  the  package  itself,  to  make  sure  the  dynamic  model  is  of 
sufficient  detail  to  accurately  desci’ibe  the  dynamic  characteristics,  or  to  point  out  the  need 
for  use  of  test  modal  data  exclusively. 

Figures  60  and  61  show  axial  and  lateral  model  predictions  for  accelei’ometei’s  SAN  and  9AN. 
It  is  obvious  that  the  correlation  is  terrible,  so  it  was  decided  to  investigate  higher  harmonic 
dynamic  models .  (Remember  that  the  axial  and  lateral  model  are  the  0™  and  1st  harmonics, 
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Figure  63.  Axial  Mode  Shape,  Test  Vs.  Analysis 

and  higher  harmonics  are  two  and  above.  This  means  that  the  circumferential  shell  dis¬ 
placement  varies  as  cos  n  9,  where  n  is  the  harmonic  number).  Ln  order  to  do  this  a  one 
inch  strip  of  the  shell  at  8AN  was  treated  as  a  ring  in  lieu  of  using  a  higher  harmonic  SABOR 
model  which  would  take  more  time.  The  frequency  associated  with  the  second  harmonic,  or 
ovaling  mode,  was  computed  at  520  Hz.  The  response  was  then  determined  and  found  to 
correspond  to  the  point  shown  on  Figure  60.  It  is  then  concluded  that  the  significant  re¬ 
sponse  peaks  on  Figures  60  and  61  are  associated  with  higher  harmonic  modes  rather  than 
the  fundamental  axial  and  lateral  models.  This  indicates  that  for  skin  or  shell  mounted 
components  higher  harmonic  models  must  be  utilized  in  addition  to  the  fundamental  modes 
to  get  accurate  response  predictions  and  resulting  skin-mounted  component  environments. 

Figure  62  shows  the  predictions  and  test  data  for  accelerometer  10AL.  This  figure  shows 
that  the  response  prediction  technique  is  indeed  accurate.  The  aft  cover  is  a  simple  struc¬ 
ture,  with  the  520  Hz  mode  being  accurately  described  by  the  dynamic  modal.  It  is  seen 
that  the  peak  response  exactly  coincides  with  the  test  response  at  this  frequency. 

Based  on  Figures  50  through  62,  it  is  concluded  that  the  response  prediction  technique  is 
accurate  and  verified.  There  are  several  suggestions  which  must  be  made  to  ensure  the 
techniques  accuracy,  A  modal  survey  (three  axes)  is  a  definite  aid  in  establishing  a  repre¬ 
sentative  dynamic  model  with  requirements  for  recording  all  accelerometers  in  each  axis  of 
excitation.  In  fact,  without  the  modal  survey  it  is  unlikely  that  the  prediction  technique 
could  be  verified  because  the  initial  dynamic  models  had  to  be  significantly  reworked  to 
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obtain  the  accuracy  required.  The  acoustic  test  pointed  out  the  need  for  including  higher 
harmonic  shell  responses  for  skin  components.  Dynamic  models  and  subsequent  response 
predictions  are  required  for  both  lateral  axes,  not  just  one  lateral  axis. 

4.  DISCUSSION 

There  are  se\eral  areas  which  should  be  elaborated  upon  as  a  result  of  the  test- analysis 
correlations.  As  pointed  out  in  the  previous  subsection,  the  developed  prediction  technique 
from  Section  IV  is  indeed  an  accurate  and  powerful  tool.  However,  it  is  based  on  finite 
clement  dynamic  models.  The  previous  oubsection  clearly  demonstrates  that  dynamic  model 
accuracy  is  of  prime  consideration  if  accurate  predictions  are  to  be  made.  The  correlation 
analysis  showed  several  shortcomings  of  the  modal  survey  in  ensuring  an  accurate  dynamic 
model.  Perhaps  the  biggest  deficiency  was  the  failure  to  monitor  the  lateral  accelerometers, 
whose  sensitive  axes  were  other  than  the  test  direction.  The  acoustic  test  clearly  showed 
axial  accelerometers  responding  at  lateral  frequencies,  which  was  not  unexpected .  How¬ 
ever,  because  of  the  failure  to  monitor  the  axial  accelerometers  during  the  lateral  modal 
survey ,  there  was  no  means  to  guarantee  that  the  dynamic  model  accurately  predicted  this 
response. 

Another  problem,  which  was  reflected  by  the  poor  correlation  above  1000  Hz,  is  that  the 
nose  force  input  to  excite  the  vehicle  resonances  did  not  excite  the  higher-order  component 
modes.  This  is  because  the  mode  shapes  associated  with  these  modes  have  little  nose  mo¬ 
tion,  thereby  resulting  in  insufficient  modal  force  in  the  survey  to  excite  them.  A  way 
around  this  would  be  to  excite  the  vehicle  at  the  locations  where  the  component  packages  tie 
into  the  shell  structure .  A  better  solution  would  be  to  run  detailed  modal  surveys  of  the 
component  packages  themselves.  This  would  guarantee  that  the  dynamic  model  is  of  suf¬ 
ficient  detail  in  the  component  areas  to  describe  their  dynamic  characteristics. 

The  correlation  analysis  of  the  previous  subsection  clearly  demonstrated  that  if  shell  ac¬ 
celeration  response  predictions  are  desired,  then  the  dynamic  models  must  include  higher 
harmonics.  The  inclusion  of  higher  harmonic  models  is  only  required  in  the  event  that  a 
small  component  is  mounted  to  the  vehicle  shell  and  the  component's  vibration  environment 
is  desired.  Both  the  acoustic  test  and  the  re-entry  predictions  demonstrated  that  the  strains 
induced  into  the  shell  structure  are  of  little  consequence  with  respect  to  shell  ultimate 
strength  capability. 

To  summarize,  the  prediction  technique  does  an  excellent  job  of  describing  re-entry  vehicle 
'•ibration  response.  However,  care  must  be  taken  to  insure  the  vehicle  dynamic  character¬ 
istics  are  accurately  porfrayed. 
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SECTION  IX 


VEHICLE  RE-ENTRY  SIMULATION  TECHNIQUES 

1.  CURRENT  AEROACOUSTIC  ENVffi03STMENT  DESCRIPTION 

The  most  current  description  of  aeroacoustie  environments  during  re-entry  of  high-beta 
vehicles  is  contained  in  Reference  23,  This  reference  documents  the  results  of  a  study  con¬ 
tract  conducted  at  GE-RESD  concurrently  with  this  study  program,  but  which  focused  on 
improving  the  description  of  aeroacoustie  loadings. 

The  acoustic  description  of  Turbulent  Boundary  Layer  (TBL)  flow  in  terms  of  Overall  Sound 
Pressure  Level  (OASPL),  Power  Spectral  Density  (PSD)  and  spatial  correlation  is  essentially 
the  same  as  that  given  in  Section  HI.  Significant  asymmetric  OASPL  distributions  about  the 
vehicle  roll  axis  do,  however,  occur  for  angles  oi  attack  in  excess  of  the  vehicle  cone  half 
angle  because  of  separated  flow. 

Transitional  flow  acoustic  environment  has  an  OASPL  of  approximately  10  dB  greater  than 
the  TBL  flow  environment,  but  retains  similar  characteristics  as  TBL  in  PSD  and  spatial 
correlation.  The  length  of  vehicle  affected  by  transitional  flow  is  approximately  equal  to  the 
local  vehicle  radius. 

Wake  flow  OASPL  is  approximately  equal  to  TBL  flow,  not  15  dB  less  than  TFL  as  defined 
in  Section  EU.  Its  PSD  is  higher  than  TBL  in  the  frequency  region  less  than  5000  IIz.  Spatial 
correlation  is  assumed  to  behave  as  previously  described. 

Data  obtained  from  recent  flights  with  severe  trajectories  indicates  that  aerodynamic  oscil¬ 
lating  shocks  can  be  set  up  on  the  vehicle  nose  if  tip-ablated  shapes  create  separated  flow. 
Published  data  (24),  (25)  indicates  that  these  oscillating  shocks  produce  very  severe  OASPL, 
30  dB  higher  than  TBL,  with  very  low  frequency  content,  less  than  100  Hz.  The  shock 
acoustic  pressures  are  fully  correlated  along  the  flow  direction  (vehicle  roll  axis),  but  fully 
correlated  over  only  one  quadrant  across  the  flow  or  around  the  vehicle  circumference. 

The  combinations  of  the  various  aeroacoustie  environments  are  depicted  in  the  schematics 
of  Figure  04  for  three  significant  flow  conditions  that  occur  during  re-entry.  During  early 
re-entry,  the  vehicle  forebody  is  completely  immersed  in  laminar  flow.  An  altitude  is  scon 
reached  where  transitional  flow  occurs  at  the  rear  of  the  vehicle.  As  the  vehicle  further 
penetrates  the  atmosphere,  the  transitional  flow  zone  moves  forward  on  the  vehicle,  with 
turbulent  flow  occurring  aft  of  the  transitional  flow  zone.  Finally,  transition  reaches  the 
nosetip  and  the  vehicle  is  completely  immersed  in  turbulent  flow.  In  addition,  at  this  point 
in  the  trajectory,  transition  causes  increased  heating  and  subsequent  ablation  of  the  nosetip. 
Depending  on  the  shape  that  the  ablated  nosetip  assumes,  separated  flow  and  oscillating 
shocks  can  occur  at  the  nose,  both  phenomena  causing  very  high  overall  sound  pressure 
levels. 
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^"isure  64,  Current  Aeroacoustic  Environment  Description 
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2.  CURRENT  .RE-ENTRY  VIBRATION  ENVIRONMENT  DESCRIPTION 

Vibration  environments  at  internal  shelf-mounted  components  measured  during  nominal  re¬ 
entry  flights  are  mild.  Figures  GO  and  GG  show  these  typically  mild  vibration  environments 
measured  during  period  of  transitional  flow  and  fully  turbulent  flow  respectively.  The  high 
levels  on  these  figures  at  very  low  frequencies  are  due  to  vehicle  rigid-body  aerodynamic 
pitching  response. 

No  data  from  vehicle  forebody  shell  and  aft  cover  (base)  are  available  for  nominal  flights. 

Data  recently  obtained,  however,  during  very  severe  re-entry  trajectories  show  extremely 
high  (4000  g's)  intermittent  acceleration  pulses  of  short  duration  (0. 1  millisecond)  occurring 
on  the  vehicle  forebody  shell  structure.  These  levels  were  estimated  by  using  triggering 
characteristics  of  high-g  switches  and  saturation  characteristics  of  vibration  sensors.  The 
cause  of  this  phenomena  is  presently  attributed,  cither  directly  (linear  response)  or  indirectly 
(nonlinear  response),  to  oscillating  aerodynamic  shocks  on  the  vehicle  ablated  nosetip.  The 
vibration  environment  occurring  during  transition  on  these  severe  trajectories  and  prior  to 
nos.etip  ablation  was  mild,  and  that  occurring  in  conjunction  with  the  intermittent  acceleration, 
pulses  could  not  be  defined  because  of  sensor  saturation. 

3.  SIMULATION  TECHNIQUES 

The  two  areas  which  must  be  investigated  in  order  to  insure  that  a  ground  test  of  a  vehicle 
will  simulate  the  re-entry  response  are  the  vehicle  dynamic  characteristics  and  the 


Figure  65.  PSD  Envelope  For  Internal  Components-Transition 
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Figure  66.  PSD  Envelope  For  Internal  Components  ~  TBL 

aeroacoustic  loading.  During  re-entry,  the  heat  shield  will  be  ablating,  the  re-entry  g 
loads  will  alter  conditions  of  preload  at  the  joint  locations,  and  the  increased  temperature 
will  reduce  stiffness  of  the  shell.  These  effects  all  combine  to  change  the  ve'’>'Te  dynamic 
characteristics  compared  to  those  of  a  ground  test  vehicle.  It  is  extremely  difficult  and 
costly  to  modify  an  unaltered  prime  vehicle  to  tine  point  where  it  has  the  exact  dynamic 
characteristics  as  a  re-entry  vehicle.  Some  things  can  be  done,  however,  to  approach  the 
ideal  condition.  If  the  effect  of  the  re-entry  g  loading  is  to  loosen  field  joints  through  reduc¬ 
ing  the  assembly  preload,  this  could  easily  be  represented  by  lower  assembly  preload. 
Although  costly,  it  is  possible  to  design  and  fabricate  a  heat  shield  which  describes  the 
predicted  re-entry  ablated  condition  of  reduced  thickness  and  modulus. 

Techniques  for  simulating  the  aeroacoustic  environments  are  shown  schematically  in 
Figure  67  for  the  three  flow  conditions  previously  described  in  Figure  64.  It  is  seen  that 
simulation  of  the  movement  of  transition  down  the  vehicle  entails  the  use  of  baffles;  this 
should  present  no  orobiem,  since  it  was  done  in  the  testing  performed  in  the  AFI'DL  Sonic 
Fatigue  Facility.  Simulation  of  the  oscillating  aerodynamic  shocks  and  separated  flow  on 
the  vehicle  •  u  requires  local  random  force  input  utilizing  a  shaker.  If  the  intermittent 
high  acceleration  pulses  described  previously  are  due  to  a  nonlinear  vehicle  response  such 
as  field  joint  or  equipment  shelf  impacting,  then  the  shaker  input  through  snug-fitting 
connectors  should  suffice.  If,  however,  the  intermittent  acceleration  pulses  are  due  to  a 
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similarly  described  aerodynamic  shock  loading  on  the  vehicle  nosetip,  then  an  impacting 
device,  such  as  loosely  fitting  connectors  at  the  nose,  must  be  used  in  conjunction  with  the 
shaker. 

Finally,  adequate  simulation  may  demand  that  the  test  vehicle  be  mounted  on  springs  such 
that  a  rigid  body  pitch  and  yaw  frequency  («20  to  40  Hz)  similar  to  re-entry  is  realized, 
since  some  components  may  be  susceptible  to  low-frequency  environments. 

4,  AFFDL  SONIC  FATIGUE  FACILITY  CAPABILITIES 

This  facility  has  several  test  chambers.  The  large  test  chamber  has  the  average  physical 
dimensions  of  42  x  50  x  70  feet,  and  can  be  operated  in  either  the  reverberant  mode  or  pro¬ 
gressive  wave  mode.  The  high  intensity  sound  in  the  test  chamber  is  generated  by  35  dis¬ 
crete  frequency  sirens  with  a  range  of  50  to  10, 000  Hz.  The  sirens  can  be  modulated  to 
produce  a  narrow  band  random  noise.  A  maximum  acoustic  power  of  one  million  watts  can 
be  radiated  from  the  main  bank  within  a  limited  frequency  range.  A  maximum  overall  sound 
pressure  level  of  174  dB  can  be  attained  near  the  sirens.  Full  scale  specimens  can  be 
tested  in  sound  fields  up  to  1G2  dB  (re  0.0002  dynes/cm^).  Seventy-two  channels  of  data 
can  be  recorded  simultaneously  by  the  facility  instrumentation  system,  with  the  capability 
of  recording  342  transducer  outputs  by  commutating  a  number  of  data  channels.  There  is 
also  a  smaller  test  chamber  scaled  to  one-fourth  the  physical  dimensions  of  the  large  cham¬ 
ber.  This  facility  is  powered  by  a  12  kilowatt  wide-band  siren  capable  of  producing  a  con¬ 
tinuous  spectrum  over  a  frequency  range  of  50  to  12,  000  Ilz.  A  maximum  overaU  sound 
pressure  level  of  160  dB  can  be  attained  at  the  siren  horn  mouth. 

As  mentioned  previously  in  various  sections  of  this  report,  the  acoustic  test  performed  with 
the  prototype  vehicle  in  the  AFFDL  Sonic  Fatigue  Facility  did  an  excellent  job  of  simulating 
responses  to  the  baseline  re-entry  environment,  even  though  the  overall  sound  pressure 
levels  were  lower  than  predicted  for  re-entry.  This  remark  can  be  made  since  the  compo¬ 
nent  responses  in  the  acoustic  test,  as  measured  by  the  accelerometers,  exceeded  the 
responses  predicted  for  the  baseline  re-entry  aeroacoustic  environment.  In  summary,  it 
can  be  stated  that  the  AFFDL  Sonic  Fatigue  Facility  adequately  simulates  re-entry,  sine  e 
the  correlation  effects  of  the  test  chamber's  semi-progressive  environment  yields  vibration 
responses  equivalent  to  re-entry,  at  overall  sound  pressure  levels  much  lower  than  re-entry 
aeroacoustics. 
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DESIGN  REQUIREMENTS,  TEST  SPECIFICATIONS 

1.  SPECIFICATION  DERIVATION  GUIDELINES 

The  basis  for  design  requirements  and  test  specifications  are  vibration  predictions  for  all 
significant  re-entry  flow  conditions  at  all  component  locations,  utilizing  the  analysis  and 
test  procedures  previously  discussed.  The  initial  task  in  specification  derivation  is  to 
overlay  the  three  axes  of  response  for  each  component,  as  one  specification  for  all  axes  of 
testing  is  desired.  If  a  component  is  used  at  more  than  one  location  in  the  vehicle,  then 
the  responses  at  all  locations  and  all  axes  must  be  overlayed.  A  test  specification  is  then 
derived  by  enveloping  the  peaks  of  the  response  overlays.  Notching  of  the  specification  is 
appropriate  only  if  the  response  overlay  peaks  occur  in  groups  at  widely  separated 
frequency  bands  as  shown  in  Figure  08.  However,  notching  is  rarely  possible  in  defining 
re-entry  vibration  specifications  because  of  the  wide  band  nature  of  the  aeroacoustic 
environment.  Rolloffs  the  test  specification  about  the  response  peaks  should  not  exceed 
12  dB/oetave  to  lie  compatible  with  test  machine  control  limitations. 

Experience  has  shown  that  enveloping  procedures  just  described,  results  in  test  overall 
levels  that  are  approximately  twice  the  response  overall  levels,  thus  providing  sufficient 
conservatism.  Additional  safety  margin  is  provided  by  the  test  duration,  which  should  be 
approximately  ten  seconds  per  axis  or  thirty  seconds  total  test  time.  Re-entry  flight  time 
from  onset  of  transition  to  impact  is  on  the  order  of  ton  seconds,  with  duration  of  highest 
level  vibration  even  less  than  that.  The  short  duration  of  test  time  will  require  the  use  of 
a  near  prime  component  simulator  to  be  used  in  shaping  the  test  specification  during 
initial  testing. 
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Figure  68.  Test  Specification  Notching 
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2.  TYPICAL  COMPONENT  TEST  SPECIFICATIONS 


Utilizing  Uk  nuincious  re-entry  response  analyses  conducted  on  CE-RESD  high-beta  re¬ 
entry  vehicle  pi'ograms  and  available  flight  data,  typical  component  test  specifications  were 
developed  using  procedures  discussed  above.  Figures  09  and  70  display  the  results  of  this 
effort  with  lest  specification  for  internal  shelf-mounted  components  and  obeli -mounted 
(forebod},  aft  cover)  components  respectively .  The  figures  provide  the  ability  to  define 
icalistic  component  re-entry  vibration  specifications  for  proposal  or  preliminary  design 
purposes,  but  obviously  more  representative  curves  should  be  derived  for  a  specific 
vehicle,  utilizing  the  analysis  and  test  procedures  outlined  in  this  report. 

An  attempt  is  made  in  Figure  71  to  give  a  best  estimate  of  acceleration  shock  pulse  environ¬ 
ments  experienced  m  the  vehicle  shell  structure  if  nose  ablation  causes  nosetip  separated 
flow  and  oscillating  shock  waves. 

3.  INTERNAL  ACOUSTIC  ENVIRONMENTS 

Previously,  in  Section  VII,  it  was  noted  that  vehicle  internal  acoustic  environments  meas¬ 
ured  in  the  testing  at  AFFDL  Sonic  Fatigue  Facility  had  overall  levels  36  dB  lower  than  the 
external  acoustic  environment.  This  means  the  external  excitation  is  fifty  times  greater 
than  the  internal  environment,  and  therefore  acoustic  excitation  of  components  in  conjunc¬ 
tion  with  vibration  testing  is  not  necessary. 
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Figure  69.  Component  Re-entry  Vibration  Test  Specification  (Internal  Self  Mounted) 
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SECTION  XI 


CONCLUSIONS  AND  RECOMMENDATIONS 


1.  VIBRATION  RESPONSE  PREDICTIONS 

The  main  objective  u£  tliia  study  contract  in  developing  a  verified  vibration  response  pre¬ 
diction  technique  has  been  achieved.  However,  the  results  of  the  test/analysis  correlate..! 
ha\e  illuminated  certain  areas  previously  not  addressed  or  little  understood.  Primarily, 
it  is  now  evident  that  aeroacoustic  environments  cause  significant  re.  _  onses  in  all  three 
axes  of  the  vehicle.  Prior  to  this  study ,  it  was  thought  that  axial  response  was  dominant 
and  the  modal  survey  and  acoustic  tests  were  planned  accordingly. 

A  very  significant  item  learned  from  the  acoustic  test  is  that  the  acoustic  environment  also 
excites  vehicle  shell  modes,  which  become  the  largest  contributor  to  shell  response. 
Finally,  test  results  reveal  that  internal  component  shelves  are  excited  at  their  higher 
resonant  frequencies  as  well  as  the  expected  lower  modes. 

Based  on  the  results  of  this  study,  the  following  is  recommended  in  predicting  vibration 
response  to  acoustic  environments: 

(1)  All  vehicle  rxes  should  be  considered  simultaneously  in  analytical  predictions. 

(2)  Shell  response  should  be  modeled. 

(3;  Modal  surveys  for  dynamic  model  verification  should  be  conducted  in  each  vehicle 
axis,  monitoring  all  axes  of  internal  responses. 

(4)  High  frequency  modeling  of  internal  equipment  shelves  must  be  included, 

(5)  Modal  surveys  of  internal  equipment  shelves  alone  should  be  done,  or  at  least 
provision  for  excitation  at  shelf  attachment  structure  during  vehicle  modal 
surveys  should  be  made. 

2.  ACOUSTIC  TESTING 

The  utilization  of  a  nose/fixture  combination  was  on  excellent  approach  for  providing  nose 
excitation  and  vehicle  suspension  during  acoustic  tests. 

It  is  certainly  evident  from  the  acoustic  test  results  that  either  type  of  chamber  configur¬ 
ation  (reverberant,  or  progressive)  is  sufficient.  This  is  because  the  proximity  of  the 
vehicle  to  the  siren  bank  precludes  good  reverberation;  hence  the  environment  tended  to  be 
progressive.  The  acoustic  environment  provided  enough  spatial  correlation  decay  around 
the  vehicle  circumference  to  cause  sufficient  lateral  response.  The  use  of  baffles  to  alter 
acoustic  environments  has  been  proven  to  be  an  effective  testing  technique. 
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In  future  acoustic  testing  it  is  recommended  that: 

(1)  more  microphones  be  placed  to  record  '’•rcumferential  variations  in  acoustic 
environment,; 

(2)  use  of  strain  gage  instrumentation  be  eliminated  because  of  extremely  low  levels; 
and 

(3)  high  and  low  frequency  portable  sirens  should  be  employed  to  raise  sound  pressure 
levels  in  those  frequency  bands. 

3.  SIMULATION  TECHNIQUES 

The  fairly  well-defined  types  of  aeroacoustic  environments  (turbulent  boundary  layer, 
transition,  wake)  are  at  best  difficult  to  simulate  in  ground  tests.  Simulation  of  oscillating 
shocks  and  separated  flow  will  require  a  better  definition  of  these  environments.  In  sum¬ 
mary,  exact  simulation  of  aeroacoustic  environments  in  ground  tests  is  virtually  impos¬ 
sible,  and  therefore  simulation  of  r<  soonses  is  most  expedient. 

The  primary  recommendation  of  this  ..  tudy  program  is  that,  now  that  feasibility  has  been 
demonstrated,  further  tests  and  analyses  should  be  performed  to  thoroughly  evaluate  the 
response  simulation  techniques  proposed  in  this  report.  This  effort  should  also  include 
tasks  to  determine  if  nonlinear  responses  in  a  vehicle  structure  could  result  in  high  inter¬ 
mittent  acceleration  pulses. 
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